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ABSTRACT 



The radiational model employed here computes the planetary 
albedo and the solar absorption by atmospheric layers and by 
earth's surface for the primitive equation model of FNWC. 
Large-scale cloud parameterization in several layers enters 
into these computations. The solar computations are made on 
a gridpoint basis employing the water-vapor mass over each 
gridpoint. Longwave cooling effects over the same layers and 
at the earth's surface are calculated using emissivity formulas 
after Sasamori, and are also dependent upon the cloud- 
parameterizations . 

Two forms of the cloud-parameterizations were tested using 
FNWC gridpoint oceanic data for 16 January 1974. The objective 
was to determine the parameterization which better verified 
the radiational balance as a function of latitude at the 
tropopause, as compared with the January satellite climatology. 
Separate heat-budget checks of the tropospheric columns and of 
the oceanic -water mass were tested for consistency with January 
climatology. The smaller cloud-parameterization values gave 
the best verification. 
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I. INTRODUCTION 



This thesis represents a continuing study of a radiative 
heating parameterization for use in appropriate o-layers of 
the Fleet Numerical Weather Central (FNWC) numerical predic- 
tion system. The immediate goal of this thesis is to analyze 
the heat balance of the earth-ocean system utilizing FNWC 
gridded data-fields at constant pressure levels for 16 
January 1974. These gridpoint data were selected along four 
meridians, three in the Pacific and one in the Atlantic 
Ocean. Most of the gridpoints were located in the Northern 
Hemisphere (Fig. 1) . 

The most important influence on the radiational disposi- 
tions (short and long-wave) of this study was the specifica- 
tion of amounts of clouds in two specific layers. The 
initial specification of the fractional amounts CL(1) and 
CL(2) are based on large-scale formulations developed by 
Smagorinsky [1960] and first used in a similar study based 
upon 16 October 1973 data by Warner [1974] . 

The governing equations of the radiational transfers of 
the model have been devised by Martin [1972, 1974], who 
modified the streams of radiation, both solar and terrestrial, 
to respond to the presence of clouds in specified amounts in 
the two layers. The radiational model developed by Martin 
and utilized here has similarities to those in use in UCLA 

i 

and NCAR General Circulation models, but with special adapta- 
tion to the o-levels of interest to FNWC. 
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Figure 1. FNY/C polar stereographic grid and meridians (lines 
1, 2, 3, and 4) selected for study. The longitudes A are 
shown for each meridian as well as the extent considered of 
each meridian. 
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Warner found that his global albedo estimates were too 
high as compared with the satellite climatology of Raschke 
et at. [1973] , thus leading to underestimates of the radia- 
tive balance both at the top of the atmosphere and at the 
surface for the data of his experiment. In this study, the 
original Smagorinsky [1960] cloud-specifications CL i and CL 2 
were computed for use in the model; however, reduced cloud- 
amounts (2/3-CL) in each layer were also included for radia- 
tive budget comparisons. 

The comparative results afforded by the "2/3-CL" model 
gave reasonably close agreement as a function of latitude 
with the radiative climatology of Raschke et at. [1973] , 
for the NIMBUS III period 21 January - 3 February 1970, 
examined by these experimenters. 

In order to analyze the effects of the "2/3-CL" versus 
the "full-CL" model used by Warner, oceanic and tropospheric- 
heat balances were also computed for both sets of cloud con- 
ditions for 16 January 1974. It was necessary in this 
connection to augment the purely radiative model by a turbu- 
lent boundary-layer model for sensible and latent heat 
transports as adapted from Kaitala [1974] . 

In the 2/3-CL case, the zonally averaged results for 
the components of the ocean- tropospheric system gave realis- 
tic results when compared to climatology for this mid-winter 
period . 

The possible error-source of using one set of computation- 
al data to be representative of mid-January climatology is 
recognized. Nevertheless, the heat-package results computed 
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with the 2/3-CL cloud model were encouraging when examined 
against climatology. The major effect noted by the reduction 
in cloud cover was the reduction in global albedo in both 
low- and mid-latitudes thus contributing to greater surface 
net heating rates. This latter result is in general accord 
with recent observations of Von der Haar and Hanson [1969], 
Von der Haar and Oort [1973] , as well as with Raschke et at. 
[1973]. However, the specification of layered amounts of 
clouds in the tropics still seems to give somewhat excessive 
global albedos. It is felt that the geometry of solar radia- 
tion streams impinging upon cumulus subgrid cells in the 
tropics must be reexamined so as to divert a greater per- 
centage of solar radiation downward, as compared with cor- 
responding cloud- layering reflective effects in middle and 
high latitudes. 
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II. DATA PREPARATION 



Temperature and humidity data utilized in this study 
were taken from gridpoint soundings along four meridians 
of the FNWC 63-by-63 gridmesh. The data gridpoints along 
these meridians (Fig. 1) were chosen so as to ensure condi- 
tions representative of oceanic locations for computations 
of the radiational and other heat budget items. The reasons 
for choosing oceanic locations for such computations were 
threefold: (1) the maritime-area heating-rate computations 

are likely to be relatively stable timewise over the data- 
day (16 January 1974); (2) the constant o = P/tt surfaces 
of the FNWC primitive equation system are close to being 
constant pressure levels; and (3) values of the heat-budget 
items computed here may be compared with corresponding cli- 
matology during mid-winter. Corresponding values over land 
areas appear to be less amenable to comparison. 

Along the four oceanic meridians used, a total of 67 
gridpoint soundings were selected from I = 1 to I = 25 on 
lines 1, 2, and 3 (over the Pacific Ocean), and 26 soundings 
from I = 63 to I = 38 along the Atlantic meridian, line 4 
(cf . , Fig. 1). The device of selecting gridpoint soundings 
chosen from the indicated meridians of the FNWC polar 
stereographic map, made it unnecessary to resort to spatial 
interpolation between original data gridpoints along the 
meridians . 

The gridpoint soundings from the data lines just mentioned 
were taken from the original FNWC 63-by-63 gridpoint analyses 
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of temperature at the surface and of T(P) and q(p) at nine 
other standard pressure levels up to and including 100 mil- 
libars (Table I (a)). Those soundings were then modeled 
into a radiative sounding model compatible in the vertical 
with FNWC primitive equation prediction levels (Fig. 2) . 

The Pacific Ocean gridpoint data cases were taken from 
the 0000GMT, 16 January, 1974 FNWC analyses, while the 
Atlantic Ocean data cases were taken from the 1200GMT, 16 
January 1974 analyses. The reason for the choice of the 
OOOOGMT time-analysis was to correspond to a set of "sound- 
ings" as close as possible to local solar noon over the mid- 
Pacific, whereas those of 1200GMT approximate a near-noon 
conditions over the Atlantic meridian. 

At each gridpoint selected, the original data listed the 
five moisture values given dewpoint depressions for the 
analysis levels from 925 mb to 400 mb. Since the standard 

instrument level vapor pressure (e - ) was missing from the 

sir 

original FNWC analysis at all gridpoints, e - was approxi- 

sir 

mated using the FNWC field of e , a computed value of vapor 
pressure at about 20 meters above MSL in the turbulent 
boundary layer. The FNWC values of T and e over the ocean 

A A 

had been made available using an operational planetary bound- 
ary layer model detailed by Kaitala [1974]. 

In Sections V and VI it was necessary to recompute values 
of e (and of T ) by application of the FNWC planetary bound- 

A A 

ary layer model to data reported in the gridpoint soundings. 

It was found that FNWC values of e^ as initially reported were 
accurate enough to be representative of the near-surface 
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Figure 2 . Five-layer radiative sounding used in this study. 
Levels are identified by their values on the k-scale, while 
layers are identified by their level boundary indices in 
parentheses, e.g. (8,10). Pressure-scaled water vapor and 
C0 2 mass increments U and C, respectively, are integrated 
with respect to the surface and introduced at even levels 
while the mixing ratio, q, is formulated at odd levels. The 
temperature, T, is expressed for all levels. Amounts of 
clouds CL(1) and CL(2) in the layers shown have been param- 
eterized for consideration of their radiative effects. 
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Table 1(a). Example of typical FNWC sounding for gridpoint 



(1,1). The humidity parameters between 925, . . . , 400mb are 
dew-point- depressions, and at the surface and top of the 



constant flux 


layer 


(9999*) are vapor 


pressures . 


Pressure 


(mb) 


T(°C) 


Humidity Parameters 


surface , 


lOOOmb 


25.8 


e =24.9 mb 

X 


925 




18.7 


1. 5°C 


850 




14.9 


1 . 6°C 


700 




5.9 


0 

n 


500 




-10.7 


4 . 9°C 


400 




-22.4 


2. 0°C 


300 




-37.3 




250 




-46.2 


- 


200 




-56.2 




150 




-67.4 




100 




-80.7 




9999* 




T =22.9 

X 





♦Code 9999 indicates data taken from the top of the constant- 
flux level (CFL) of the FNWC initial data program. 
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Table 1(b) . Example of the corresponding radiative sounding 



with mixing ratio listed at odd k-levels (Fig. 2). Addition 
ally, water vapor and C0 2 absorber masses, cloud amounts, 
CL(1) and CL(2) where applicable in the radiative sounding. 



Pressure (mb) 


T(°K) 


Mixing Ratio 
g/kg 


1000 


25.8 


15.49 


900 


17.1 


12.76 


800 


12.1 




700 


5.9 


6.12 


600 


- 1.3 




500 


-10.7 


2.31 


400 


-22.4 


1.37 


300 


-37.3 


.54 


200 


-56.2 




100 


-80.7 


.02 


0 


-80.7 


.00 



Absorber Masses Cloud Amounts 



Water Vapor 
(gm/cm 2 ) 


C0 2 

(cm/ cm 2 ) 


CL(1) ,CL(2) 


2.392 


45.53 


CL (2) 
1.000 


3.349 


83.53 


CL(1) 

.650 


3.633 


113.35 




3.679 


134.00 




3.680 


138.67 




3.680 


143.35 
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humidity values, e . , even though the computations for e 

made in Sections V and VI employed other parameter- values 
for these computations. 

In addition to the parameters listed in Table I (a, b) , 
surface geostrophic windspeeds were also calculated. In the 
procedure for calculation of the windspeeds at the gridpoints, 
use was made of the 100 mb FNWC analysis to determine the 
contour heights as D- values. At certain locations, climatic 
geostrophic windspeeds were used (Atkinson, 1970). A 
more complete explanation of the use of the geostrophic 
windspeed calculations and its use in connection with the 
planetary boundary layer sensible and latent heat transports 
is found in Section V, and the windspeed values are listed 
in that section. 

To perform the radiative calculations indicated in Table 
I (b) , data-processing methods were applied to convert both 
the temperature and mixing-ratio data at k-levels to obtain 
the cloud amounts, CL (1) and CL (2) at certain required k- 
level boundaries (Fig. £) . It should be observed that all 
soundings in this study start at sea level with the approxi- 
mation of surface pressure ir = 1000 mb. Thus the k-levels 
correspond approximately to the FNWC a-levels and are taken 
as P = 1000., 900., 800., ...., 200., 100., 0.0 respectively. 

A. INTERPOLATIVE PROCESSING TO k-LEVELS IN RADIATIVE SOUNDINGS 

1 . Temperatures 

At each mandatory level of Table 1(a) between the 
1000.,..., 100 mb, gridpoint temperatures were listed. An 
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assumption was made that at the top of the atmosphere, the 
temperature was isothermal from 100 mb to p = 0.0 mb. The 
listed FNWC sea-surface temperature was taken as T.^. At 
the remaining k-levels, air temperatures were derived from 
either the corresponding FNWC listed temperature or where 
necessary by interpolation from listed values of the original 
FNWC data. The interpolation scheme used was that of the 
three-point Lagrangian formula in the vertical. 

T k - T o ((P k' P l )(P k' P 2 ))/t(P o" P 1 ) (P 0 ' P 1 ;i) 

* T i aPk-vfVppvctPj-vcPi-pp) (2-D 

+ T 2 ( ( P k -P 0 )(P k -P 1 ))/((P 2 -P 0 )(P 2 -P 1 )). 

It is noted that in Equation 2-1, k is the required pressure 

level needed for the T, interpolation. The level P lies 

k r o 

above and the pressure levels P^, and P£ lie below level 
k. The temperatures T q , T-, , and T ^ are at the same levels 
as P Q , P^ and respectively. 

2 . Humidities At k-Levels 

Moisture parameters from the original gridpoint 
soundings (Table I (a)) were converted into mixing ratios 
at each of the original sounding levels (Table I (a)). 

Then utilizing Equation 2-2 the near-surface vapor pressure 
(taken as e ) was employed to calculate the mixing ratio at 
k = 10 as follows: 

q 1000 = 621.97 (e air /1000) (2-2) 

where e a ^ r = e x is in mb and Q^qqq is in gm(kgm) . The 
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remaining mixing ratios were calculated from the original 
FNWC dewpoint depression data levels (Table I (a)) using 
Equation 2-3 (after Fleagle and Businger, 1963): 

q(P) - « - t-9t«a» -. t*--B/Tl exp [ , (2-3) 

where q(P) = mixing ratio at pressure level P and 
A = 21.656 
B = 5418. 0°K 

T = temperature (deg. K) at level P 
Tp = dewpoint (deg. K) at level P 
T-Tp = dewpoint depression at level P. 

The computed q-values were subsequently interpolated to k- 
levels using the procedure indicated in Equation 2-1. An 
example of the interpolated q-values is shown for the case 
of the radiative sounding at point (1, 1) (Table I (b)). 

3 . Humidity Extrapolation to p < 300 mb 

Due to the fact that most radiosonde humidity data 
for p < 300 mb is either unreliable or not available, a 
power-law extrapolation formula (Equation 2-4) similar to 
one used by Smith [1966] , was employed to obtain q-values at 
k = 3, 2, and 1: 



qfP) = 

<*5 




(2-4) 



where q(P) = mixing ratio at level k in g Kg 

q^ = mixing ratio at level P = 500 mb in g Kg 1 
P = pressure at level k in mb. 

Equation 2-4 was solved for the best-fit parameter relating 
the variables y = log (q (P) /q^ ) and x 5 log(P/500). Six q- 
values from the original sounding including P = 1000, 925, 
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850, 700, 500, and 400 mb from Table I (a) are used in the 
best fit solution (e.g., Crow et al. t I960): 



X 



if 6 

i = l 



y^ x 



i i 



i = 6 

,E 1 x. : 

1=1 l 



(2-5) 



After the determination of A from Equation 2-5, extrapolation 
of q-values from P = 400 to P = 300, 200, and 100 mb, was 
performed utilizing Equation 2-4. 

The procedure of the previous paragraph was carried 
out at each gridpoint. As a test of the usefulness of the 
X-profile technique in determining q-values above 300 mb, a 
correlation coefficient was statistically computed for each 
X-profile. The correlation coefficient computed was between 
y and x and had the form: 




(2-6) 



As was found in the studies of Jenks [1974] and Warner [1974] , 
the correlation coefficient was found generally to lie in 
the range of .95 to .99. Such high correlations indicate 
that the X-profile technique for determining upper atmosphere 
q-values has realistic estimation value for the vertical-scale 
involved. 

4 . Pressure-Scaled Absorber Masses 

By using the mixing-ratio values computed at the 
five odd k-levels, the layer pressure-scaled water vapor 
absorber masses denoted by AU (Fig. 1) were calculated for 
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each layer centered at the five odd k-levels. The general 
equation for computing the pressure-scaled water vapor mass 
was : 

^k 72 

U(k-l,k+l) - — t— [ 13YJ725 ]' 72 (2-7) 

where the pressure ratio raised to the exponent 0.72 has 
been employed (after Moller and Raschke, 1964). k is the 
odd-level centered in each vertical layer concerned and 
AP = 200 mb was used for each of the five layers. The pro- 
cedure was followed for each layer (8, 10), (6, 8), (4, 6), 

(2, 4), and (0, 2). 

The algorithm for computing the integrated water 
vapor mass above the earth's surface was then developed as 
follows : 



U(10) 


= 0.0 






U (8) 


= U (10) 


+ U (8 


, 10) 


U (6) 


= U(8) 


+ U (6 , 


8) 


U (4) 


= U(6) 


+ U (4 , 


6) 


U(2) 


= U(4) 


+ U (2 , 


4) 


U(l) 


= U(2) 


+ 0.5U 


(0, 2) 


U(0) 


= U(2) 


+ U(0, 


2). 



In the computation of the carbon-dioxide scaled 
absorber mass, a computational routine similar to the one 
used for water vapor was employed. Because it is customary 
to state the CO ^ absorber masses in terms of the N. T. P. 
volume of the gas in cms over a cm 2 of the earth's surface, 
the CO 2 absorber mass (AC) centered on an odd k- level becomes 
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(2-9) 



AC k = C(k-l,k+l) - 3.14 x 10“ ( |£- )( lQ ^ 25 )°' 72 

Here k is the odd-level indicator, 3.14 x 10 4 is the con- 
stant mixing ratio of CC >2 by volume, while AP/gp s is the N. 

T. P. thickness of any layer of pressure- increment AP = 200 
mb. P s is the constant density in a homogeneous atmosphere 
of depth H, in which 



T(Surf) = 273.16 deg. K and = 1013.25 mb. 

Finally, Equation 2-9 is shown to be represented in the form 



C (k-1 ,k+l) 



3.14 x 10 (H) ( 



AP 

1015.25 



)( 



1013.25 



)' 



72 



RT 

H = — - = 7.9954 x 10 5 cm. 
g 



(2-10) 



Equation 2-10 was then utilized for each odd level k = 9, 7, 
5, 3, 1. 

As was done in the water- vapor case, an algorithm 
similar to Equation 2-8 is developed to find the integrated 
carbon dioxide mass above the earth's surface: 



C(10) 


= 0.0 






C (8) 


= C(10) 


* C(8, 


10) 


C(6) 


= C (8) + 


C (6 , 


8) 


C (4) 


= C(6) + 


C (4 , 


6) 


C (2) 


= C(4) + 


C (2 , 


4) 


C(l) 


= C(2) + 


0 . 5C (0 , 2 


C(0) 


= C(2) + 


C(0, 


2). 



( 2 - 11 ) 
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B. CLOUD PARAMETERIZATION 



The relative humidities and thus the saturation vapor 
pressure at levels k = 5 and k = 9 are used in the calcula- 
tions of the fractional cloud cover amounts in layers (4, 6) 
and (8, 9). The saturation vapor pressures at each level 
were computed using the following equations (cf.. Equation 



2-3) : 

621.97 exp [ A - ] 

q s (5) 500 5 — t 2 ' 12 

621.97 exp [ A - y— ] 

q s (9) * goo — ( 2-12 (b)) 

where A = 21.656 and B = 5418.0 deg. K. 

The relative humidities were then calculated from: 



R.H.(S) = (2-13 (a) ) 

R.H.(9) = 3^- . (2-13 (b)) 

Using these relative humidity values, the fractional cloud 
amounts for each layer were parameterized employing equations 
similar to Smagorinsky [1960] : 

CL ( 1) = 2.0 (RH (5) ) - 0.7 (2-14 (a)) 

CL (2 ) = 3.33 (RH (9) ) - 2.0. (2-14 (b)) 

Equations 2-14 have been employed in the radiative studies of 

Jenks [1974] for a spring day in 1973 and by Warner [1974] for 
a case in mid-October, 1973. In the current study, Equations 
2-14 are tested in a winter sample of FNWC data for 16 
January 1974. 
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Smagorinsky ' s parameterization of CL(1) and CL(2) per- 
mits cloudiness fractions of 1.30 and 1.333 respectively with 
RH = 1.0. Values of CL of this magnitude were allowed orig- 
inally by Smagorinsky to suggest the type of supersaturation 
which accompanies precipitation. However, in all of the 
radiative models in use (that of FNWC, UCLA, NCAR, etc.), 
it was considered correct procedure to limit CL < 1.0 in 
both layers. In the radiative calculations of this study an 
arbitrary reduction by one third of that computed by the 
Smagorinsky formulations (2-14) was made so that neither 
fractional cloud-cover could exceed unity. For radiative 
calculations, the degree of supersaturation in clouds was 
considered not of importance. 

Thus CL(1) and CL(2) computed by (2-l4a, b) are 
replaced respectively by their comparison cloud-covers, 



CL (1) ' and CL(2)', defined by 

CL (1) ' = 2/3 CL (1) 
CL (2) ' = 2/3 CL (2) 



(2-15) 



computed at each gridpoint. The objective for testing the 
reduced cloud-covers (2-15) was to determine the relative 
sense of the displacement of the radiative balance of the 
earth atmosphere system. Such radiative tests were especially 
crucial since evidence from Warner's [1974] tests made with 
unmodified (CL(1), CL(2)) indicated that the planetary albedo 
was larger than observations from NIMBUS III (Raschke, et al. 
1973). The one-third reduction in the Smagorinsky calculated 
amounts used in this study is an initial effort to tune cloud 
amounts to give radiational results in closer agreement with 
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the latest satellite results such as are found in Raschke 
et at. [1973]. Also, the fractional cloud amounts were 
considered functions of large-scale effects only. Small 
scale convective activity, seasonal conditions (other than 
the soundings themselves) , and latitudinal effects were not 
considered in specifying this reduction factor. The purpose 
of the reduced cloud-parameterization modification was for 
estimating radiational effects only. 



C. WINDSPEED USED IN TURBULENT HEAT TRANSFERS 

For all gridpoints located between 5°N and 5°S, gradient 
windspeeds were taken from the climatological data for the 
month of January (Atkinson, 1970). Gradient isotach values 
were utilized at the tropical gridpoints noted due to the 
small Coriolis parameter f in Equation (2 - 16) (below) in the 
latitudes between 5°N and 5°S. 

At all other gridpoints, Equation (2-16) was utilized to 
compute geostrophic windspeed given by 



IV 



g 1 



= £ 
f 



JAlL 
2 A 

v nr 



(2-16) 



here g = 9.80665m sec 



m = map factor = 



1.86603 



<t> 

d 



1 + S in 

latitude 

381. x 10 3 m = nominal FNWC grid-dimension 



f = ( 



4tt 

86,400 



) S incf) 



(2-17) 



|Az| = /(Zi, - Z 2 ) Z + ( 2 3 - zTT 2 . 



(2-18) 
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In (2-16) and (2-18), Az is the contour- di fference cen- 
tered in general on the gridpoint. At the map edge grid- 
points of meridians 1, 2, 3, 4 it was possible to compute 
only a forward differenced version of Equation 2-18 (cf . , 
Fig. 3). Thus at edge-points "0", Equation 2-19 was used 
(below). Note, the grid distance d/m in Equation 2-19 is 
one-half of that used in Equation (2-16). 




g 

f 



h z l 

(d7i) 



(2-19) 



with 

Az = /(zi» - z o ) 2 + ["zl - z o ) 2 . (2-20) 



All geostrophic windspeeds were then converted to surface 
windspeeds using the following empiricism due to Langlois 
and Kwok, [1969] . 

V s = .8 | V g | + 2.2 (2-21) 

where = surface windspeed (mps) 

V = 1000 mb geostrophic windspeed. 

g 

Note that in the interval 5 to 10 degrees of latitude in 
either hemisphere, the value of the Coriolis parameter was 
arbitrarily set at f = .25 x 10 **sec 1 , following Kaitala 
[1974], while for | <p | > 10°, the actual value of f for this 
particular latitude was used. 
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Figure 3. Section of FNWC polar stereographic meridian il- 
lustrating the method of obtaining contour gradients in the 
vicinity of a gridpoint. The values of contour-heights 
correspond to the 1000 mb D-values in the neighborhood of 
the gridpoint. 
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III. TERRESTRIAL RADIATION 



A. THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas for computing emissivities for ter- 
restrial flux calculations in the atmosphere were developed 
by Sasamori [1968] in association with the operation of the 
NCAR General Circulation Model. These empirical formulas 
for water vapor and CC^ emissivities were matched to the 
theoretical values used in the Radiation Chart developed by 
Yamamoto [1952]. This chart was proved by Jenks [1974] to 
be accurate enough to simulate the empirical methods of 
Sasmori and is used here as a schematic guide for the inte- 
gration of the radiative transfers through the various layers 
in the radiative soundings described in Section II. 

For use in the FNWC heating package, the essential long- 
wave flux formulas required are the following: 

F 1q * = net IR flux at earth, k = 10 

Fg* = net IR flux at level, k = 6 

F 2 * = net IR flux at level, k = 2 

F610 = net IR flux divergence in the layer (6, 10) 

F26 = net IR flux divergence in the layer (2, 6). 

In order to get the IR flux-divergences (i.e., cooling 
rates) in the layers (6, 10) and (2, 6), the differences 
F^* - F^O* an< ^ ^2* " ^6* mus ^ computed. The detailed 
scheme for making such computations for various combinations 
of cloud cover CL(1), CL(2) is similar to that used by Warner 
[1974] and is explained in the following subsections B, C, D, 
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and E along with the diagrams Figures (4a,b,c,d), (5a,b,c,d) 
and (6a,b,c,d). A complete listing of net IR fluxes and 
cloud amounts was computed at each gridpoint. An example 
of this listing is shown for "full-cloud" amounts, Equation 
2-14, in Table II. One must further have a physically sound 
representation of the emissivity e as a function of both 
water-vapor and CO 2 absorber masses along the sounding, that 
is, along the path of integration in order to make IR net- 
flux calculations. For further explanation of the emissivity 
formulas (after Sasamori 1968) used in the quadrature scheme, 
refer to Appendix A. 



B. NET FLUX F 1Q * 

As depicted in Table I (b) , radiative soundings have been 
computed as the combination of parameters U(k, 10), C(k, 10), 
T^, for each required level k. Also, listed at each grid- 
point are the cloud-cover parameters CL(1) and CL(2) com- 

4 

puted from (2-14). The grid area may then be thought of 
being composed of areal fractions (or weights) W(0, 0), 

W(l, 1), W(l, 0), and W(0, 1) which are defined by Equations 
3-1, j3-4 as follows: 

(a) totally clear skies, case (0, 0) 

W(0, 0) = (l-CL(l)) (1-CL(2)) (3-1) 

(b) overcast in both layers, case (1, 1) 

W(l, 1) = (CL (1) ) (CL (2) ) (3-2) 

(c) overcast in upper layer only, (case 1, 0) 

W(l, 0) = (CL (1) ) (1-CL (2) ) (3-3) 

(d) overcast in lower layer only, (case 0, 1) 

W(0, 1) = (1 -CL (1) ) (CL (2) ) . 
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CL( 1 ) = 0.0 CL( 2 ) = 0.0 WT = 0.0 F 10 * = .1639 F 6 * = .2701 F 2 * = .3575 ly/min. 
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By multiplying the appropriate weighting factors by 



the respective net flux computations F^q*( 0, 0),..., 

Fio*(0, 1) the composite net fluxes at level k = 10 for 
the combined cloud (CL(1), CL(2)) coverage over the grid 
area are then obtained. Such composite net flux values 
have the following form using level k = 10 as an example: 

F 1q * = W(0, 0)F 1 q *( 0, 0) + W(l, 1)F 1() *(1, 1) 

+ W(l, 0)F 1 q *(1, 0) + W (0 , 1)F 1(J *C0, 1). 

1. Net Flux with Clear Skies, F ]_q *(0, 0) 

• In Fig. 4 (a) the hatched area represents downward 

flux F104- at k = 10 or in integral form: 

A 10 „ 

F wc* - J 0 E »c dB ■ B K ■ ° T K 

B^ is the Stephen-Boltzman blackbody flux at T^. e is the 
combined water-vapor and CO^ emissivity along the sounding, 
and is related to absorber-masses as detailed in Appendix A. 

The unhatched area therefore depicts the net flux at k = 10 
with' clear skies given in the equation 

yOlO 

F 10*C°’ °5 - A U • W dB > B 10 = St * T io'- ( 3 ‘ 6) 

B- U 

2. Net Flux with Overcast Clouds in Both Layers, F ^q *( 1,1) 
Figure 4 (b) illustrates this case. By using the 

trapezoidal summation techniques, F104- at k = 10 is depicted 
once again by the hatched area representative of the case of 
an overcast with base at k = 9. The unhatched area depicts 
the integral 



(3-7) 



F * f 1 
10 u » 



1 ) 




.e ) dB 
wc ; 



with overcast clouds at both upper and lower levels (Fig. 4 

(b)). 

3. Net Flux with Overcast Clouds in Upper Layer Only, 

F 10 *(l. °) 

Here in Fig. 4 (c) , the hatched area depicts the 
total F104- emitted down from the cloud base at k = 6 , plus 
an increment in the upper right corner of Fig. 4 (c) , which 
represents the contribution to F4- from the water vapor and 
CO 2 in the atmospheric layer from k = 6 to k = 10. The un- 
hatched area depicts the integral 

p 10 

F 10*d- °> - A (1 ' £ wc^ dB - C 3 - 8 > 

-D “ U 

4. Net Flux with Overcast in Lower Layer Only F ^q *( 0, 1) 
This case is the same as depicted for case F^q*(1, 1) 

in Fig. 4 (d) and is explained in Section III.(B).2. e wc is 

the temperature - dependent emissivity for the completion of 

the integration from level k = 1 to the origin (B = 0 , e = 

1.0) of the Yamamoto Radiation Chart, applied for T < 210 

deg. K. e wc is the temperature - independent emissivity applied 

where T > 210 deg. K in the level from k = 1 to k = 10. For 

the specific formulations of c and e after Sasamori [1968] 

^ wc wc 

see Appendix A. 

5. The Composite F ^q * (CL(1) CL(2)) Equation 

The principle involved in compositing net fluxes has 
already been expressed by Equation 3-5. The following equa- 
tion can be thought of as the general form of the quadrature 
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scheme used to evaluate the integrals for computing F^q* for 
any combination of (CL(1), CL(2)). It is also used in com- 
puting the integral for the composite cloud-cover case F^q*, 
corresponding to (3-5): 

F 1q *(CL 1 ,CL 2 ) = [l-CL(2)]{(B 10 -B 6 )-.5[e wc (8,10)(B 10 -B 8 ) 

+ C^ c (8,10) +£wc (6,10))(B 8 -B 6 )]} 

+ (1-CL(2)) (1-CL(1)0{B 6 -.5[(£ wc (6,1O) 

+ £ wc (4 ’ 10))(B 6- B 4 ) + ^wc^ 4 ’ 10 ) C 3 ' 9 ) 

+ e wc (2,10))(B 4 -B 2 )+(e wc (2,10) 

+ e wc (l,0))(B 2 -B 1 ) + £ wc ((0,10),T 1 )*B 1 ]} 

+ CL(2){(B 10 -B g ) [1-.5 £ wc (9,10)]}. 

By setting CL(1) = CL(2) = 0.0, Equation 3-9 reduces 
to the form needed in calculating F 4 q*( 0, 0). Similarly, one 
can obtain the equation for computing F^q*( 1, 1) as depicted 
on the chart in Fig. 4 (b) by setting CL(1) = CL(2) = 1.0. 
Finally, the other two subcases F^q*( 1, 0) and F 4 q*( 0, 1)) 
can be formulated by setting CL(1) = 1.0, CL(2) = 0.0, and 
CL(1) = 0.0, CL(2) = 1.0, respectively in Equation 3-9. 

Similar procedures were followed for computing net IR flux 
at levels k = 6 and k = 2, as described in the Subsections 
C and D (Figs. 5 (a), (b) , (c) , (d) and 6 (a), (b) , (c) , (d) , 

respectively . 



C. NET FLUX F 6 * 

1 . The Composite F ^. * Equation 

The following expression, which represents the F^* 
composite net flux, may be thought of as the weighted sum of 
the unhatched areas of Figs. 5 (a), (b) , (c) , and (d) , with 

the appropriate weight factors W(0, 0), W(l, 1), W(l, 0), 
and W(0, 1) of Equations 3-1, 2, 3, and 4 applied. The re- 
sulting expression for Fg*(CL^, CL^) becomes (after Warner, 

1974) 

F 6 * = [l-CL(l)]{B 8 -.5[e wc (6,8)(B 8 -B 6 )'+E wc (4,6)(B 6 -B 4 ) 

+ (^ wc (4,6) + e wc (2,6))(B 4 -B 2 ) + (c wc (2,6) + c wc (1,6))(B 2 -B 1 ) 
+ £ wc ((0,6),T 1 )*B 1 ]}+(1-CL(1))(1-CL(2)){(B 10 -B 8 ) [1. 

- .5(e wc (6,8) + c wc (6,10))]}+CL(l){(B 8 -B 6 )* (3-10) 

[l-.5e wc (6,8)]} + CL(l) (1-CL (2) ) { (B^-Bg) * 
[l-.5(e wc (6,8)+e wc (6,10))]}. 

2 . Net Flux With Clear Skies F^*(0, 0) 

As was described in Section III.B.5 concerning the 
use of the composite case (Equation 3-9), Equation 3-10 can 
be reduced to give expressions for F^* for the various cloud- 
cover cases (0, 0), (1, 1), (1, 0), and (0, 1). The resulting 
expressions are depicted as the unhatched areas in Figs. 5 (a), 
(b) , (c) , and (d) . For the clear sky case of F^*, Equation 
5-10 would be employed with CL(1) = CL(2) = 0.0. This is 
depicted by the unhatched area in Fig. 5 (a). 



3 . Net Flux With Overcast Clouds in Both Layers F^*(l, 1) 
For this case, Equation 3-10 would be used with CL(1) = 

CL(2) = 1.0. This is depicted by the unhatched area in Fig. 

5 (b). 

4. Net Flux With Overcast Clouds in Upper Layer Only 

yq, Q) 

The expression needed in this case, and which is 
depicted by the unhatched area in Fig. 5 (c) , is formulated 
by using Equation 3-10 and setting CL(1) = 1.0 and CL(2) = 

0 . 0 . 

5. Net Flux With Overcast Clouds in Lower Layer Only 

y to, i) 

The expression needed in this case, representable by 
the unhatched area in Fig. 5 ( d) , is formulated by using 
Equation 3-10 and setting CL(1) = 0.0 and CL(2) = 1.0. 

D. NET FLUX F 2 * 

A procedure analogous to the ones used in the computa- 
tions of F.. * and F,-* net fluxes can be used in calculation 
10 o 

the net flux at level k = 2. F ? * is represented pictorially 

by the weighted sum of the unhatched areas of Figs. 6 (a), 

(b) , (c) , and (d) . Once these areas have been computed with 

the proper weighting factors applied, as was done similarly 
in Equations 3-9 and 3-10, the following general expression 
for the composite result for F 2 * results 

F 2 * = [l-CL(l)]{B 8 -.5[e wc (2,4)CB 4 -B 2 )+(e wc (2,4) 

+ £ wc( 2 ’ 6 )HB 6 -B 4 ) + (e wc (2,6) + c wc (2,8))* 

(B 8 -B 6 )+e wc (l,2)(B 2 -B 1 )+e wc ((0,2),T 1 )*B 1 )]} (3-11) 
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(3-11) 



+ (1-CL(1))(1-CL(2))*{(B 10 -B 8 ) [ 1- . 5 ( £ wc (2,8) 



+ e wc (2,10))]} + CL(l){B 4 -.5[£ wc (2,4)(B 4 -B 2 ) 



+ e wc (l,2)(B 2 -B 1 )+e wc ((0,2),T 1 )*B 1 ]} 



Equation 3-11 can also be derived in the computation of the 
net fluxes for the general cloud-cover cases once the proper 
weight factors for the simplified net fluxes are introduced. 
These simplified F 2 * cases correspond to the (0, 0), (1, 1) , 

(1, 0), (0, 1) cloud-overcast/clear combinations and are 
represented by the notation F 2 *(0, 0), F 2 *(l, 1), F 2 *(l, 0) , 
F 2 *(0, 1), respectively. The corresponding chart depictions 
are displayed in Fig. 6 (a), (b) , (c) , (d) , respectively. 

E. APPLICATION TO HEAT BALANCE COMPUTATIONS 

In order to compute the heat balance of the atmosphere, 
layers (2, 6) and (6, 10) are treated as subject to solar 
heating (Section IV) . In order to complete the atmospheric 
radiational balance, the long-wave mean cooling rate in the 
FNWC model was computed as follows (after Kaitala [1974]): 



where the combination 26 and 610 indicates the layer boundaries 
Equations 3-12(a) and (b) represent the cooling effects by ter- 
restrial radiation in the indicated layers, tvhich is applicable 
at the time of the analysis. The radiative balance at the sur- 
face is computed using the absorbed solar insolation (Section 
IV.B-C) and the net flux F^q* from Equation (3-9). A complete 
heat balance at the surface is computed by subtraction of the 



F26 = F 2 * - F 6 * 



(3-12 (a) ) 



F610 = F * - F * 
roiu r 6 r 10 



(3-12 (b)) 
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evaporative and sensible heat transport (Section V) from the 
surface radiative "balance." All these balances are dis- 
played in cross sectional display forms in Section VI. The 
IR-cooling analysis was carried out at each gridpoint and 
the cooling rates were stored by gridpoint along each meridian. 
The zonally-averaged atmospheric balances are detailed in 
Sections VII and VIII for 16 January 1974. 

F. STATISTICAL RESULTS AND COMPARISONS 

A statistical test of the F^g* numerical results (Equa- 
tion 3-9) was conducted by utilizing a linear regression 
program for regressing the 93 values of F^g* on the cor- 
responding sample values of B^g and B-^g/e as predictors. The 
program used was based on the BMD02R in the Biomedical set 
of programs (after DIXON, 1973) . The resulting form was 



where B in = Stephen-Boltzman blackbody radiation at the sur- 
face with temperature T^g. B^g = St * T^q 4 



a and b in Equation 3-13 are the desired regression coeffi- 
cients and are dependent on B^ and B-g/cT (after Brunt, 1932). 
The predictand used for the regression was F^q*, which in- 
cluded the downward contribution of CO 2 as well as that of 
water vapor in the computation of net flux at level 10. The 
statistical equation found in this study is Equation 3-14 
where R m is the multiple correlation coefficient. 




(3-13) 



e = surface layer vapor pressure (mb) . 



F 1q * = .3564 B 1q - .0224B 10 /e 

R m = .9955. 
m 



(3-14) 
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The high value of the R shows that B iri and B iri /e are good 

m 10 10 

predictors o f F 10 * and therefore of obtained from 
F 1q * = B 1q - F^. Introduction of F-^q* into 3-14 permits 
solution of the latter for the downward flux at the surface. 

F^ is given by 

F d = B 10 ( ( 1 - a ) -b/e} 

or (3-15) 

F d = B io^' 6436 + • 0224/e) . 

Equation 3-15 is to be compared with a similar result after 
Swinbank [1963] 

F d = B i 0 C* 64 + - 037v/ ®)- (3-16) 

In summary, the high value of the correlation coefficient of 
Equation 3-14 along with a strong similarity between Equations 
3-15 and 3-16 seems to indicate that Equation 3-14 may be 
regarded as an oceanic version of the Brunt downward- f lux 
equation. Also, it is noted that F^ depends most strongly 
upon B 1q and only slightly on (B^g/e"). This was to be ex- 
pected over the oceanic regime where nearly constant relative 
humidities exist over large, expansive areas. 

A second statistical test was based upon a possible re- 
lationship between cloudy sky and clear sky cases of F^g*. 

In this test, the ratio F^* (CL) /F-^q* (0 , 0) was used as the 
predictand. The predictor was an approximation to the com- 
posite total opaque cloud at each gridpoint, defined as: 

CL = CL (1) + CL (2) ( 1 - CL (1) ) . (3-17) 

The regression was run twice, first utilizing Fjg*(CL) com- 
puted using the amounts of CL(1) and CL(2) as calculated from 
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the Smagorinsky formulas (Equations 2-14 (a), (b)) and second, 
utilizing F^q*(CL) computed from the two- thirds amounts of 
CL(1) and CL(2) . The resulting regression formulas based 
upon a best-fit to the 93 soundings used in this study were 
listed in (3-18, 3-19) 



FULL CL: 


F 10 *(CL) 


= F 1q *(0, 0) (1- 


. 75CL) , R = .9916 


(3-18) 


2/3 CL: 


F 1q *(CL) 


= F 1q *(0, 0) (1- 


. 7 3CL) , R = .9936. 


(3-19) 



A comparison equation after Sverdrup [1942] is included as 



(3-20) 



F 1 q *(CL) » F 1q *(0, 0) (1 - . 83CL) . (3-20) 

The high correlation coefficients found with Equation 3-18 and 
3-19 along with the similarity to the Equation 3-20 (after 
Sverdrup, 1942), indicates the general capability of the IR- 
radiative model employed in this study to account for the 
effective net radiation at the surface, F^q* (CL^ , CL^) , both 
in the completely clear and cloudy-sky cases. 

The mean values of the parameters in Equations 3-18 and 
3-19 over the 93 gridpoint soundings were as follows: 

FULL-CLOUD STATISTICS REDUCED-CLOUD STATISTICS 

F 1q *(CL) » .0749 iy min" 1 F 1Q *(CL) = .0950 SLy min' 1 

F 10 *(0, 0) = .1518 ly min' 1 E 1Q *(0, 0) = .1518 ly min" 1 



CL = .6951 



CL = .5190 



In the above table, superior-bar symbols denote sample means. 
These results indicate that the surface net flux is decreased 
by the respective ratios ( . 0 769/ . 15 18) = .507 with CL = 

.6951 and (. 0568/ . 1518) = .3740 with CL = .5190. These results 



53 



indicate that the downward flux F^gl(CL) is, on the average, 
50.7 percent larger than the clear-sky case utilizing the 
cloud-parameterization model (Equations 2-14(a), (b)) in 
this study with full cloud amounts. On the other F^qI(CL) 
for the two-thirds cloud-model, the downward flux is increased 
by only 37.4 percent relative to the clear-sky case. 

G. COMPARISONS WITH CLIMATOLOGICAL RESULTS OF F 2 

Comparison was made of F 2 by the model-computations with 
satellite measurements of total long-wave flux to space after 
Raschke et al. [1973] for the NIMBUS III period 21 January - 
3 February 1970. F ? is used instead of F 2 * since the satel- 
lite does not measure the downward flux at k = 2 . F z is 
computed using Equation 3-11, which was also used for F 2 * 
but omitting those emissivity- terms involving levels above 
k = 2, that is e (2, 1) and e (2, 0), from Equation 3-11. 
These terms are associated with downward flux at k = 2. The 
results for 16 January 1974 were computed both for the full- 
cloud and the two-thirds cloud-amount cases, based upon all 
gridpoint soundings. These F 2 - values were then interpolated 
to whole multiple of 5° latitude along each of the four merid- 
ians, and the results latitudinally averaged to get a "zonally- 
averaged" value valid for each 5° latitude band used in this 
study . 

Table III (a) shows the zonally- averaged F 2 - values com- 
pared with those after Raschke et al . For consistency in 
comparison of the F 2 results deduced here with those reported 
by Raschke, the latter's results were obtained by interpola- 
tion from charts using the same oceanic meridians as those 
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Radiative Statistics 



Cloud Statistics 



F 2 F 2 F 2 FULL-Cloud 2/3-Cloud 

Lat (FULL-CL) (2/3-CL) (Raschke) CL AVG. CL AVG. 



20S 


.311 


.325 


. 370 


. 932 


.758 


15 


.321 


. 334 


.378 


. 927 


.733 


10 


. 344 


. 354 


.377 


.867 


.634 


5 


.354 


.363 


. 347 


.746 


. 536 


0 


. 361 


.369 


.351 


.790 


.556 


5 


.371 


.378 


.368 


.766 


.527 


10 


.382 


. 387 


. 388 


. 557 


.376 


15 


. 386 


.393 


. 394 


. 532 


.342 


20 


.390 


. 388 


.391 


.443 


.279 


25 


.361 


. 362 


. 371 


.463 


. 320 


30 


. 354 


.328 


. 344 


. 559 


.411 


35 


. 336 


. 342 


. 319 


. 516 


.382 


40 


. 315 


.319 


.304 


.628 


. 529 


45 


.278 


.270 


.292 


.833 


.740 


50 


.264 


.280 


. 270 


.911 


. 759 


55 


.217 


. 233 


. 255 


. 931 


. 749 


60 


. 209 


. 226 


.258 


.764 


. 528 


65N 


. 235 


. 235 


. 230 


0.0 


0.0 


Wt. 

Avg. 


. 333 


. 338 


. 344 


.688 


. 514 


Table 


III (a) . 


Comparison 


of zonally- 


averaged 


longwave flux 


to space, F 2j 


as found by 


this study 


for both 


full-CL and 


2/3-CL 


cloud 


cases for 16 


January 1974, and by 


Raschke et al. 



(1973) based upon NIMBUS III measurements. Also included are 
composite cloud amount fractions used in the two cloud 
models of this study. Flux values in ly min -1 . 
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used in this study. For further insight into the compari- 
sons, the cloud statistics, CL and 2/3-CL values are listed 
at each latitude, where CL is given by Equation 3-17. In 
the bottom line of each column in Table III (a) is listed 
the cosine -weighted mean of each set of the column values 
over the listed latitude range. The resulting oceanic-means 
of F 2 obtained by the various techniques are as follows: 

F 2 (Full cloud amount this study) = .333 Zy min 1 

cloud amount this study) = .338 ly min 1 

F 2 (after Raschke) = .344 Jly min 1 . 

Note that the Raschke results do not correspond to speci- 

fic known CL - values. The only conclusion that can be drawn 
from Table III (a) is that both cloud-amount systems give F 2 
results which are reasonably close to the values reported by 
Raschke. This is especially true in the subtropics, 15°-20°N, 
where peak values of F 0 are found, and in high latitudes where 
minimum values of F 2 (after Raschke, 1973), are found. By 
these results alone, neither the full-cloud model nor the 
2/3 CL can be selected as preferable. The limitations of 
the comparisons made here are obvious, when it is realized 
that at latitudes between 20°S-5°S and between 60°N-65°N in 
Table III (a), fewer than 4 meridional lines enter into the 
zonal averages listed. For all other latitudinal average 
values, four meridional lines were used in this averaging. 
These same limitations will apply to Table III (b) below. 

Note finally that the model results F 2 (full) and F 2 (2/3 CL) 
are based upon a single-day's observations, 16 January 1974. 



56 



Lat 



F i o * 

(FULL CL) 



F i o * 
(2/3 CL) 



F i o * 
(Budyko ) 



20S 


. 054 


.078 


.091 


15 


.050 


.079 


. 092 


10 


. 055 


. 086 


.083 


5 


.063 


. 090 


. 085 


0 


.052 


.082 


.084 


5 


.056 


. 086 


. 080 


10 


. 085 


.107 


. 086 


15 


.094 


.115 


.094 


20 


. 102 


.120 


.102 


25 


.104 


.120 


.098 


30 


. 096 


.111 


.099 


35 


.107 


.122 


.102 


40 


. 096 


.108 


.109 


45 


.056 


. 067 


.098 


50 


.054 


. 069 


.099 


55 


.052 


.068 


.102 


60 


.085 


.106 


.104 


65N 


. 095 


. 095 


.095 


Wt. Avg. 


.075 


.096 


.094 


Table 111(b). 


Comparison of 


zonal ly- aver aged 


net fluxes at 



the surface F i0 * as found by this study for both full-CL and 
2/3-CL cloud cases for 16 January 1974, and comparative 
values after annual climatology of Budyko (1956). Flux 
values in ly min“ 1 . 
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H. COMPARISONS WITH CLIMATOLOGICAL RESULTS OF ¥ 1Q * 

Table III (b) depicts the zonally-averaged values of 
the F 1q * at the earth’s surface obtained in this study using 
both full-cloud model and the two-thirds cloud model. A 
comparison of the model -computed F^^-values for the full 
and two-thirds cloud cases is then made with annual mean 
values of obtained from climatological studies (Budyko, 

1956) . The F^q* values shown in Table III (b) were computed 
from Equation 3-9 using the full CL^,CL 2 values or the re- 
duced (2/3CL^, 2 / 3 CL 2 ) in the respective cloud-models. 

Since Budyko' s mean values of F-^o* were based on annual cli- 
matology, his overall latitudinal range of F^q* values is 
not as great as that of F^q* by either cloud-model presented 
in Table III (b) for 16 January 1974. 

It is to be noted that the locations of highest values 
of F 1q * after Budyko agree with those computed for both F^q* 
(CL) and F^q*( 2/3 CL), namely in the subtropics (latitude 
20°N) , and again in the latitude range 35-40°N. The "hemi- 
spheric" weighted mean of F^q* = .096 ly min 1 by the re- 
duced-cloud model is slightly greater than the cosine-weighted 
mean of Budyko (.094 ly min *). However, the weighted mean 
of F^q*(CL) resulting from the full-cloud model is consider- 
ably smaller, F^q* (CL) = .075 ly min 1 . The chief differences 
in F 10 *(CL) as a function of latitude compared with that of 
Budyko is accounted for by the marked minimum in F^q* cor- 
responding to the relatively high composite CL-amounts in 
low-latitudes (from 20S to 5N) and again in the latitude band 
45-55°N where in both areas typical minimum values of F^q*(CL) 
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= .055 Zy min 1 occur. Such a minimum of i- n the regions 

45-55°N is generally associated with the location of the 
polar-front zone. However, in this same region the 2/3-CL 
model gives average values of F^q* = .068 Zy min . The 
zonal-mean values of CL and 2/3-CL applicable to Table III 
(b) by latitude are exactly as listed in Table III (a). It 
is not clear from the comparisons of Table III (b) whether 
the higher values of F^q* deduced from the 2/3-CL model 
represents a closer approximation to reality than the F^q*(CL) 
values resulting from Equation 3-9 with the full-cloud model 
values . 

It is to be expected that analogous comparisons made 
using the solar-insolation model (summarized in Table IV) 
will indicate some relative advantage for using the full or 
two-thirds cloud model. 
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IV. SOLAR RADIATION ■ 



A. PARTITION OF SOLAR INSOLATION 

The solar constant assumed in this study at level k = 0 
(top of atmosphere) was 2.00 ly min 1 (Joseph, 1971). This 
flux-value was further depleted by 4% to account for the 
attenuation caused by oxygen and ozone above the tropopause. 
This left the value S = 1.92 SLy min 1 at level k = 2 to be 
used in this study as the effective solar constant . 

Equation (4-1) was then utilized to compute the effective 
solar insolation at the tropopause (k = 2) as follows: 

F(2) = S [ ]" 2 Cos z (4-1) 

m 

where S = effective solar constant at k = 2 

Cos z = cosine of the zenith angle for the Julian date used 

r/ r = ratio of the actual earth-sun distance to the mean 
in 

earth-sun distance for the Julian date used in 
this study. 

The Smithsonian Meteorological Tables (List, 1958) gives the 
ratio r/r^ and the solar declination 6 for 16 January, 0000GMT, 
as these values are employed. These values are, respectively: 

— = .983715 
m 

6 = -21.08 degrees of lat. 

6 is used in evaluating the cosine of the solar zenith angle, 
given by 

Cos z = Sin <f> Sin 6 + Cos <J> Cos 6 Cos h (4-2) 
where <J> is the latitude and h is the hour angle of the sun 
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relative to the meridional data lines. For example, Fig. 1 
makes it clear that for 



line 1 h = 55° 

line 2 h = 10° 

lines 3, 4 h = 35° 

at the times of the synoptic charts (0000GMT and 1200GMT, 
respectively) . Sin (j> was computed according to the standard 
polar stereographic projection formula applicable to the 
FNWC base chart given by 



r 2 - [ (I ~ 32) 2 + [ (J- 32) 2 ] 

Sin <f> = — (4-3) 

r E 2 + [ (I -32) 2 + [ CJ-32) 2 ) 

where r^ 2 = 973.752. Thus Sin <t> assumes the following func- 
tional form in terms of the FNWC grid-coordinate I (Fig. 1) : 
Lines 1, 3, 4 



c . . 973. 752-2(32-1) 2 

<p 973.752+2(32-1)2 



(4-4 (a)) 



Line 2 



C- A 973.752- (32-1) 2 

bin 9 973. 752+(32-I) 2 



(4-4 (b)) 



or, conversely I is given in terms of 4> by 

Lines 1, 3, 4 I = 32 - 22.065 [ 1 (4-4 (c)) 

Line 2 I = 32 - 31 . 205 [ — (4-4 (d)) 

1,..., 25 for Lines 1, 2 
I = 9,..., 25 for Line 3 

63,..., 38 for Line 4. 

For lines 1, 2, and 3 the gridpoint soundings correspond to 
0000GMT 16 January when solar noon occurs at the 180^ 
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meridian. Line 4 gridpoint soundings correspond to 1200 GMT, 
16 January, when solar noon was at the Greenwich meridian. 

A very simple partition of solar insolation was utilized 
in this study after Joseph [1971] . It consisted of dividing 
the insolation F(2) into two parts at level k = 2. One part 
was considered to include all wavelengths A > .9 pm where 
absorption by water vapor and carbon dioxide bands are the 
most prevalent attenuation processes in clear air. This part 
of solar insolation was listed as F (A) energy and considered 
subject to water-vapor absorption but not to Rayleigh scat- 
tering. For those wavelengths A < .9 pm, absorption of the 
solar insolation energy by water vapor was considered neg- 
ligible. This part of the solar insolation band was listed 
as F(S) and was subject only to Rayleigh scattering attenua- 
tion in clear air. The two partitions are formulated as 
follows : 



In this study, the introduction of two cloud decks produced 
cloud-reflectivity effects upon both the F (A) and F(S) solar 
energy insolations. However, in the clear areas around any 
gridpoint the absorption-attenuation only applies to the 
F(A) insolation while the Rayleigh scattering-attenuation 
applies to the F(S) insolation only. 

B. DISPOSITION OF F(S) INSOLATION 

In the treatment of the F(S) insolation, Joseph [1971] 
determined that Rayleigh scattering reflectance in clear skies 



F (A) = . 349 F (2) 

F (S) = .651 F (2) . 



(4-5 (a)) 
(4-5 (b)) 
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(after Coulson, 1959) could be effectively approximated by 
least squares in the following form 



a (R) = .085 + .25074 [log ( J- Sec z )] (4-6) 

o 

where P q = 1013.25 mb. In Equation 4-6, tt/P o = 1 in view of 
the fact that mean sea level pressure tt is close to 1000 mb. 
Also 

Sec z = (Cos z) 1 



with Cos z given by Equation 4-2. 

The surface albedo a(G) is another reflective parameter 
utilized in this study. Over oceanic areas the following 
formula for a(G) after Gates et al. [1971], was utilized: 

a(G) = max (.06, .06 + .54 (.7 - Cos z)}. (4-7) 

1 . Clear Sky Case 

In the clear sky (0, 0) case the F(S) insolation was 
subjected only to Rayleigh scattering reflectance a(R) and 
to surface reflectance a(G) . Considering the likelihood of 
a succession of multiple reflections between earth and 
atmosphere, the F(S) insolation actually penetrating the 
earth's surface after scattering is given by 



IS10(0,0) = F(S) [l-a(R)] [l + a(R)a(G) + . • • (a(R)a(G)) n 

+ • • -]*(l-a(G)) (4-8 (a)) 



that is, by 

IS10 (0 , 0) = F(S) [l-a(R)] [l-a(G)]/[l-a(R)a(G) (4-8 (b)) 

2 . Cloudy-Sky Cases 

In the three cases in which clouds were present, F(S) 
insolation absorbed by the ground at each gridpoint was 
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computed using a variation of the following equation (after 
Arakawa, 1972) : 

IS10 (1 , 1) = F(S) (l-R(l)) (1-R(2)) (l-ct(G)) 

*{l- [(R(l)R(2)+R(2)a(G)+R(l)a(G) 

+ 2R(l)R(2)a(G)] } _1 . (4-9) 

As indicated by the notation (1, 1), denoting CL(1) = CL(2) = 

1.0. Equation 4-9 is the formula used in calculating F(S) 
insolation absorbed by the ground in the case where overcast 
clouds are present at both levels of Fig. 2. Also in Equa- 
tion 4-9, constant cloud- ref lectance values were chosen, 
namely R(l) = .54 for the mid-level clouds between k = 4 
and 6, and R(2) = .66 for the low-level clouds between k = 

8 and 9. Both cloud-reflectance values are as suggested by 
C. D. Rodgers [1967] . 

For the other cloud cases, the following changes 
were applied to Equation 4-9. In the (1, 0) case (CL(1) = 

1.0, CL(2) = 0.0), consider that R(2) = 0 in (4-9), from 
which it follows that 

IS10(1, 0) = F(S) (l-R(l)) (l-a(G))/[l-R(l)a(G))] . 

(4-10) 

In the case (0, 1), R(l) = 0.0 in (4-9) so that 

IS10 (0 , 1) = F(S)(l-R(2))(l-a(G))/[l-R(2)a(G)]. 

3 . Composite F(S) Insolation 

Equations (4-8), (4-9), (4-10), and (4-11) were 
utilized in the computation of the cloud-weighted F(S) in- 
solation penetrating the earth's surface considering the 
areal -weights of the cloud combinations denoted by (0, 0), 

(1, 1), (1, 0) and (0, 1) about a gridpoint. The resultant 



F(S) insolation penetrating the earth's surface denoted by 
IS10 is therefore expressible as 



IS10 (CL (1) , CL (2) ) = IS10(0, 0) W(0, 0) 



+ IS10(1, 1) W(l, 1) 

+ IS10(1, 0) W(l, 0) 

+ IS10(0, 1) W(0, 1) . 



(4-12) 



Here the weighting factors W(0, 0), W(l, 1), 1V(1, 0) and 
W(0, 1) are as computed in Equations 3-1, 3-2, 3-3, and 3-4 
respectively. Note finally that the part of F(S) insolation 
reflected to space is found by subtracting IS10(CL(1), CL(2)) 
from F (S) . 

C. DISPOSITION OF F (A) INSOLATION 

The fractional portion of the solar insolation subject 
to absorption by atmospheric water-vapor and carbon dioxide 
are covered in the following subsections. 

1 . The Clear-Sky Case (0, 0) 

The Manabe-Moller absorptivity function provided the 
necessary absorptivity values for the key layers in this case. 
The form of this absorptivity function is 



Here a_ is the absorptivity applied to the pressure-scaled 
water vapor mass between levels 2 and k (Fig. 2) along the 
zenith slant-path angle z. The resultant absorbed insola- 
tional energy in the particular layer (2, k) is then given 
by the Manabe-Moller relation 



a(2, k) = . 2 7 1 [U (2 , k)Sec z]’ 303 



(4-13) 



A (2 , k) = 0 . 2 71F (A) [U (2 , k) Sec z]’ 303 



(4-14) 
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The two layers of interest in which absorption was computed 
were (2, 6) and (2, 10). The absorbed insolation in the 
layer (6,10) was then computed by 

A (6 , 10) = A (2 , 10) - A (2 , 6). (4-15) 

It is noted here that the water-vapor mass above level 2 was 
assumed negligible in this study. 

By subtracting A(2, 10) from F (A) , the direct transmis- 
sion of F(A) insolation impinging at the earth's surface was 
determined. The transmission of F (A) insolation is then fur- 
ther reduced by the transmissivity (l-a(G)), after surface- 
reflectance which leads to the earth-absorbed result 

IA10 (0 , 0) = F(A) (1- . 271[U(2, 10)Sec z] * 303 } (1-a (G) ) . 

(4-16 (a)) 

The transmitted energy impinging upon the earth iust prior to 
absorption is 

TRANA (0 , 0) = IA10 (0 , 0) / [1-«(G)]. (4-16 (b)) 

2 . Cloudy Cases 

In order to compute meaningful dispositions of F(A) 
insolation in cloudy-sky cases, cloud reflectivities and 
cloud absorptivities after C. D. Rodgers [1967] were utilized. 
The ref lectance -values used here are different from those 
suggested by Rodgers for the F(S) wavelengths. The cloud 
reflectivities used here are RA(1) = .46 and RA(2) = .50. 

In this case there are also cloud-absorptivit ies to be con- 
sidered. These were taken as A(l) = .20 and A(2) = .30, 
respectively. In the following discussions the cloud con- 
ditions are considered totally overcast in each of layers 
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indicated when the notations ((1, 1), (0, 1), and (1, 0)) 

are utilized. 

Schematic representations of the computations per- 
formed in the various cases ((1, 1), (0, 1), (1, 0)) are 

displayed in Figures 7, 8, and 9. These figures suggest the 
parameters required in computing the disposition of incoming 
F(A) - insolation from level k = 2 to the earth's surface 
(k = 10) . 

Figure 7 depicts the disposition of F (A) insolation 
in the overcast case with both high and low cloud layers 
Cl, 1). Figure 8 shows the case (1, 0) with - an upper over- 
cast only. Figure 9 displays the case (0, 1) for a low 
overcast. The equations (after Warner, 1974) which relate 
to the parameters in these figures are listed in Appendix 
B-l, B-2, B-3, respectively. 

In Figs. 7, 8, and 9, the insolations, A24, A46, 

A68, A89, and A910, etc., represent the insolation absorbed 
in the layers concerned. Symbols F2, F4, F6, F8, and F9, 
etci, depict the streams of insolation passing through the 
indicated level. A vertical arrow implies the direction of 
insolation passage, i.e., 4 denotes downward insolation, i 
upward-reflected insolation, and 4-4- downward-reflected inso- 
lation. The absorption quantity A(6, 8)14-, for example, in- 
dicates absorbed energy remaining in (6, 8) from a downward 
reflected beam. 

Since multiple reflections occur between the earth's 
surface and a cloud base, or between cloud-layers in this 
model, the insolational amount reaching the lowermost reflecting 
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Figure 7 Schematic representation of F(A) insolation 
disposition in the case of two overcast layers. 
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Figure 8 Schematic representation at F(A) insolation dis- 
position with an upper overcast layer only. 
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Figure 9 Schematic representation of F(A) insolation 
disposition with a lower overcast layer only. 
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surface is affected by the downward reflected beams. Addi- 
tional insolational amounts remaining after more than two 
reflections are negligible and are not included in the models 
(Appendix B) . Also note that insolation reflected upward 
from a lower interface, either cloud-top or ground, to the 
base of a higher cloud deck was not subjected to absorption 
by or transmission through the cloud. This simplification 
resulted in slightly reduced F(A) insolation- ref lectance to 
space . 

A final point to note is that there are three con- 
tributions to the absorption of solar insolation within a 
layer between any two reflecting surfaces. In the Figs. 7, 

8, 9, these contributions are identified by arrows which in- 
dicate the portion of path being crossed. 

For each of the overcast-cloud combinations, the 
transmitted F(A) insolation arriving at the earth's surface 
may be defined using the following notation (TRANA) as 



TRANA (1, 1) 


= F101 


+ F1011 






(4-17) 


TRANA (1 , 0) 


= F101 


+ FlOll 






(4-18) 


TRANA (0, 1) 


= FI 01 


+ F1011 






(4-19) 


where the right side 


parameters of 


(4-17), 


(4-18), 


and (4-19) 


are derived in Appendices 


(B-l), (B 


- 2) , and 


(B r 3 ) . 




The F (A) insolation absorbed by the 


earth 


in each 


case is derived by 












IA10 (1 , 1) 


= F101 


(l-a(G) ) 


+ FI 0 1 1 




(4-20) 


IA10 (1 , 0) 


= F101 


(1 -a (G) ) 


+ F1011 




(4-21) 


IA10 (0 , 1) 


= F101 


(1 -a (G) ) 


+ F1011 


• 


(4-22) 
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In the last three equations, the quantity F10LL is small 
enough in each case, that no further reflections from the 
earth were considered. 

3 . Composite F (A) Layer-Absorptions and Surface - 

Absorption Insolation 

As has been previously discussed, the standard grid- 
area weighting scheme of this study was applied to obtain 
composite values of the absorbed F (A) insolation in key 
layers and also within the earth's surface. The weighting 
factors applied to the corresponding overcast -combination 
absorption quantities provided the following composite 
results : 

A26 (CL (1) , CL (2) ) = A26 (0> Q) W(0, 0) + A26 (lj 0) W(1, 0) 

+ A26 (0j 1} W(0, 1) (4-23) 

+ A26 (l, i)W(l, 1) 

A6 10 (CL ( 1) , CL (2) ) = A610 t0 Q) W(0, 0) 

+ A610 (1> 0) W(1, 0) (4-24) 

+ A610^ 1) W(1, 1) + A610 (q 1) W(0, 1) 

IA10 (CL (1) , CL (2) ) = IA10 (0 , 0)W(0, 0) 

+ IA10 (1 , 0)W(1, 0) 

(4-25) 

+ I A10 (0 , 1)W(0, 1) 

+ IA10 (1, 1)W(1, 1) . 

The weighting factors W(0, 0),...,W(1, 1) were first listed 
in Equations 3-1, 3-2, 3-3, 3-4. 

4 . Absorptivity (ABA) by Layers 

Here the (fractional) absorptivity as well as the 
actual insolation values absorbed in the layers are considered. 
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In the computation of absorptivity, the total undepleted in- 
solation at the top (k = 0) is used. The following equation 
was utilized in this calculation: 



- 2 



FADJ = 2 . 00 (r/r ) Cos z. 
^ ' nr 



(4-26) 



ABA, or the absorptivity of the troposphere was computed from 
the ratio of the insolation absorbed in the troposphere to 
the insolation incident at the top of the atmosphere rather 
than at k = 2: 



ada = A (2 , 6) + A (6 , 10) 

A A FADJ 



(4-27) 



D. ALBEDO (ALB) OF THE EARTH-TROPOSPHERE SYSTEM 

In considering the planetary albedo, the absorptions of 
the earth- troposphere system in both the F(A) and F(S) in- 
solational regions must be recalled by the program. Thus 
we define the reflected insolational energy in F(A) as REFA, 
and it is computed by 



REFA = F (A) - A26 - A610 - IA10 . 



(4-28) 



Likewise REFS, the reflected part of F(S) is defined by 



REFS = F (S) - IS10. 



(4-29) 



It is understood that REFA, REFS are computed for appropriate 
CL^, CL 2 at each gridpoint sounding considered. 

Summing REFA and REFS gives REF: 



REF = REFS + REFA. 



(4-30) 



Finally the planetary albedo is related to FADJ through 



ALB = 



REF 

FADJ 



(4-31) 
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E. COMPOSITE ABSORPTIVITY (ABG) BY THE EARTH - SURFACE ; 

COMPOSITE ATMOSPHERIC TRANSMISSIVITY (ATRAN) 

1 . Absorptivity (ABG) of Earth 

By summing the weighted values of F(S) and F (A) por- 
tions of the incoming insolation entering the earth, the 
total insolation absorbed at the earth's surface was com- 
puted. This quantity when divided by the extraterrestrial 
insolation gave the fractional absorptivity (ABG) of the 
earth's surface. The equation for ABG was 



ABG = 



I A10 + IS10 
FADJ 



(4-32) 



where IA10, IS10, and FADJ were defined previously by Equa- 
tions 4-25, 4-12, and 4-26, respectively. 

2 . Transmissivity (ATRAN) of the Troposphere 

Also computed was the total insolational energy 
TRAN , incident at the earth's surface just before absorption 
by the surface. This calculation is given by 



TRAN = TRANA + [IS10 / (l-a(G)]. (4-33) 

Here a(G) is as previously defined in Equation 4-7. TRANA 
is as defined by the weighted value of TRANA(0, 0), TRANA(1, 1) , 
TRANA(1, 0) and TRANA(0, 1) given by Equations 4-16, 4-17, 

4-18, and 4-19. Also note that the four terms of IS10 of 
(4-8), (4-9), (4-10), and (4-11) have the common factor 
(l-a(G)) in the numerator and therefore each F(S) insolation 
component at the earth just before absorption need only be 
divided by (l-a(G)). TRAN may thus be viewed as the total 
insolational energy incident at a pyrheliometer located at 
earth. The (fractional) transmissivity of the troposphere 



74 



(ATRAN) is then computed from 



ATRAN = TRAN / FADJ. (4-34) 

Note finally that the major dispositions of the total in- 
solation at the indicated map times have been identified by 
the fractional values, ALB, ABA, ABG, and ATRAN, each of 
which is a fractional value representing the reflectivity 
(albedo), absorptivity, or atmospheric transmissivity as 
the case may be. 

3 . Computational Check 

The computational scheme utilized in this model was 
checked by summing the fractional values ALB, ABA, and ABG 
relative to the troposphere at each gridpoint. The value in 
each case must total .96, since as previously noted the at- 
tenuation of solar insolation was taken as four percent as 
it streamed through the stratosphere. 



F. STATISTICAL ANALYSIS 

In order to substantiate some of the computations per- 
formed in this section, several of the most important items 
computed were statistically analyzed using linear regression 
computer programs from the BMD set of statistical programs 
(Dixon, 1973) . 

1 . Clear Sky Cases 

Using ALB (0, 0, z) , which is the clear sky case of 
albedo, as the predictand, and log-^Sec z and its square as 
the predictor, the following best-fit equation resulted 



ALB (0 , 0, z) = .05705 + .54529 log 1Q Sec z 

- . 18885 (log Sec z) 2 . 



(4-35) 
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.9918 



where represents the multiple correlation coefficient. 

This result was as expected over the ocean where 
a(G) and a(R) involved logarithmic dependence on Sec z. 
Average values of ALB and log^Sec z were 

ALB (0,0) = .20350 log 1Q Sec z = .34049. 

Other clear sky regression tests made use of the water-vapor 
mass path length (M) defined by 

M = (U Sec z)' 5 log 10 (U Sec z)^. (4-36) 



This parameterization of water vapor mass is similar to that 
developed by Hanson [1971] , who used a similar M as a pre- 
dictor in his empirical formulations of ABA for both clear 
and partly cloudy sky cases. In this study the two regres- 
sions attempted using M as a predictor were with ABA(0, 0, M) 
and ATRAN(0, 0, M) as predictands. The best-fit equations 



resulting were 



ABA(0, 0, M) 
R 

m 

ATRAN (0 , 0, M) 
R 



m 



.11478 + .04012 M 
.9675 

.77794 - .06337 M 
.8575. 



(4-37) 

(4-38) 



The means of ABA(0, 0, M) and ATRAN (0, 0, M) for the original 
93 soundings were 



ABA(0 , 0, M) = .15154 



ATRAN (0, 0, M) 
M 



.71987 

1.83259 (gm cm 2 )^. 



76 



2 . Relationships Between Albedo, Atmospheric Absorptivity , 
At mospheric Transmissivity, and Ground Insolation in 
the Cloudy and Clear Sky Cases 

In addition to the clear sky value of ALB, it was 
also possible to compute the composite cloudy sky case. A 
regression was formed showing the relationship between the 
ratio ALB (CL (I) , CL (2) ) /ALB (0 , 0) as the predictand and the 
total opaque cloud cover (CL) as the predictor. CL in this 
model is specified by Equation 4-39 below: 

CL = CL (1) + CL (2 ) - CL(1)CL(2). (4-39) 

CL, it was felt, gave a good approximation of the effective 
cloud-cover by the two layers of cloud amounts CL^ and 
by Equation 2-14 used in this model. As explained earlier 
in this study, a one-third reduction of both CL(1) and CL(2) 
after initial determination by Equation 2-14 was also tried 
in all computations involving cloud amounts. This test 
represented an initial attempt at tuning the cloud model for 
radiative calculations as will later be substantiated, and 
gave results in closer agreement with the latest literature 
on satellite reflectances. Thus all subsequent regressions 
involving CL will also include those best-fit equations for- 
mulated using a CL based also upon a one-third reduction of 
CL(1) and CL(2). A small letter subscript 'a' in the equa- 
tion number will indicate those computations using full-cloud 
amounts while a small letter subscript ' b ’ will indicate those 
utilizing a one-third reduction in the cloud model amounts. 

The first regression equations tested were those for 
ALB(CL) and ABA(CL) for the respective cases of "full" and 
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two-thirds cloud cover. The results are given in (4-40), 

4-41) . The symbol R m once again signifies the multiple 
correlation of the statistical regressions. 

ALB (CL (1) CL (2) = ALB(0,0)[1 + 3.59683CL 

cl 

- 1.7006CL 2 ] (4-40 (a)) 

R = .8755 
m 

ALB b (CL(l) ,CL(2) = ALB (0 , 0)[1 + 4.13667CL 

- 3. 18972CL 2 ] (4-40 (b)) 

R = .8876. 
m 

The sample means of the values in Equations 4-40a and 4-40b 
are 



Equation 4-40a 




Equation 4-40b 




ALB (CL (1) , CL (2) ) = 


. 4298 


ALB(CLCl) , CL (2) ) = 


.3640 


ALB (0 , 0) 


.2035 


ALB (0 , 0) 


.2035 


CL 


.6951 


CL 


.5176 



The correlation coefficient in Equation 4-40 appears to show 
strong dependence of ALB upon CL. Also note that ALB(CL(1), 
CL(2)) values computed here are slightly higher than reported 
by Raschke et at ., [1973], (ALB = .30) for essentially the 

same set of gridpoints during the NIMBUS III period 21 January 
through 3 February 1970. However, the one-third reduction-CL 
appears to give much closer albedo results, suggesting that 
the full cloud amounts of Smagorinsky by 2-14 are too high. 

A second regression computed was between the ratio of 
tropospheric absorptivity of the cloudy air to that for the 
clear sky case, versus CL. The following best-fit equation 
resulted: 
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ABA (CL (1) ,CL(2)) = ABA(0, 0)[1 + .5043CL 

3 . 

- . 1326CL 2 ] 

R = .8989 
m 



ABA b (CL(l) ,CL(2)) = ABA(0, 0)[1 + .5522CL 

- . 24 51CL 2 ] 

R m = .8873 
m 

Mean statistics for this regression case are 



(4-41 (a)) 



(4-41 (b)) 



Equation 4-41(a) 




Equation 4-41(b) 




ABA a (CL(l),CL(2)) = 


.1916 


ABA b (CL(l) ,CL(2) = 


. 1810 


ABA(0 , 0) 


. 1515 


ABA(0 , 0) 


.1515 


CL 

a 


.6951 


CL b 


.5176 



In (4-41), it is seen that the model specifies an increase in 
atmospheric (solar) absorptivity with increasing the cloud 
cover. This result is in agreement with studies of Plante 
[1972], Warner [1974], and with the earlier study of London 
[1957] . 

An analogous statistical regression was then developed 
for the cloudy-sky transmissivity relative to the clear-sky 
transmissivity. The resulting regressions were 



ATRAN (CL(1) , CL (2) ) = AT RAN (0,0) [1-.6020CL 

3 - 

- . 0285CL 2 ] 

R = .9895 
m 



(4-42 (a)) 



ATRAN b (CL(l) ,CL(2)) =ATRAN(0,0) [1 

+ . 1246CL 2 ] 

R = .9862. 
m 



6833CL 

(4-42 (b)) 
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Equation 4-42b 



Equation 4-42a 
ATRAN CL(1) , CL (2) = .4440 ATRAN, CL (l) , CL (2) = .4957 

3. D 

ATRAN(0, 0) = .7199 ATRAN (0 , 0) = .7199 

If Equation 4-42b is terminated at the first power of CL, the 
result reduces to 

ATRAN (CL (1) , CL (2)) = ATRAN (0, 0) [1- . 5961CL] (4-43) 

R = .9855. 
m 

The result (4-43) is in reasonable agreement with the Savino- 
Angstrom result which is 

ATRAN (CL(1) ,CL(2)) = ATRAN (0 , 0) [ 1 - (1- k) CL] (4-44) 

Here k is a slowly increasing function of latitude, ranging 
from 0.35 to 0.50, according to Budyko [1956]. 

3 . Some Conclusions 

The statistical analyses performed seem to show agree 
ment witli observational results of other investigators. Thus 
the radiation model, with the solar seasonal effects imposed 
appears to be realistic. Further, the one-third reduction in 
amounts given by the cloud parameterization model appear to 
give even better agreement with the latest satellite findings 

G. ALBEDO COMPARISONS WITH PUBLISHED RESULTS 

The tropospheric albedo computations of the solar-inso- 
lation model of this section have been presented for both 
the full-cloud and two-thirds CL cases, respectively. These 
albedos have been interpolated to whole multiples of 5° 
latitude between 20°S and 65°N. These computations were made 
for each of the four oceanic meridians and the resulting 
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zonally-averaged albedos are presented as a function of 
latitude in Table IV for both the full-CL and two-thirds CL 
cases, where comparison is also made with satellite clima- 
tology of Raschke et at. [1973] , pertaining to the NIMBUS 
III averaged albedo measurements over the period 21 January - 
3 February 1970. For sake of consistency the Raschke albedo 
analyses have been derived by interpolation to latitude using 
the same four meridians before applying the zonal - averaging 
process . 

The model results under the headings "FULL-CL" and "2/3- 
CL" correspond to the composite cloud amounts previously 
listed in Table 11(a) according to the equation 

CL = CL 1 + CL 2 - CL 1 * CL 2 (3-17) 

with either (CL^, CL 2 ) by Equation 2-14 or (CL^, CLp by 
Equation 2-15. 

The most unusual results to be observed in Table IV are 
the small albedo-values (~ .2) between 20°S to 20°N after 
Raschke, as contrasted with the corresponding values (in the 
range .3 to .4) by 2/3-CL cloud model , and more markedly with 
those of the FULL-CL model (values in the range .35 to .5). 

The "hemispheric" cosine weighted mean albedos are as follows: 

FULL-CL, ALB = .46 

2/3-CL, ALB = .38 

RASCHKE -ALB = . 30 

At every latitude (<J> < 50°) under comparison the 2/3-CL 
albedo computation is closer to Raschke 's observations than 
is the full-CL value. At latitudes <f> > 55°, there is evidence 
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ALBEDO ALBEDO ALBEDO FULL-CL (2/3-CL) 

Lat (FULL-CL) (2/3-CL) (Raschke) AMOUNT AMOUNT 



20S 




.5462 


.4458 


.2195 


.932 


.758 


15 




.5380 


.4336 


. 2050 


.927 


.733 


10 




.4901 


.3859 


. 1875 


.867 


.634 


5 




.4334 


. 3432 


. 1967 


.746 


. 536 


0 




.4447 


.3478 


.2325 


.790 


.556 


5 




.4290 


.3362 


.2225 


.766 


. 527 


10 




. 3532 


. 2859 


.2050 


.557 


.376 


15 




. 3517 


.2895 


.2075 


.532 


.342 


20 




. 3390 


. 2868 


.2450 


.443 


.279 


25 




.3664 


.3151 


.2925 


.463 


. 320 


30 




. 4217 


.3668 


.2725 


. 559 


.411 


35 




.4145 


. 3703 


.3375 


.516 


. 382 


40 




.4595 


. 4275 


. 3925 


.628 


. 529 


45 




. 5558 


. 5151 


.4250 


.833 


.740 


50 




. 5682 


. 5186 


. 5325 


.911 


.759 


55 




.5717 


.5218 


.5929 


.931 


.749 


60 




.5180 


.4848 


.5600 


.764 


. 528 


65N 




.4294 


. 4153 


.7000 


0.0 


0.0 


Wt. 

Avg 




.4622 


.3820 


.3010 


. 688 


.514 


Table 


IV 


Comparison of planetary 


albedo as found by this 


study 


for both full 


and two- 


thirds 


cloud models 


for 16 


January 


1974 and by 


Raschke 


et al . 


(1973) based 


upon NIMBUS 



III measurements. Also included are the globally weighted 
mean values at the base of the columns and the composite 
cloud amounts for both cloud cases. 
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of added surface-ice reflections contributing to the Raschke 
values of planetary albedo. This is especially true at the 
latitude 65°N where Raschke's results indicate ALB(65°N) = 

.7, whereas the model -computations with essentially zero 
cloud cover is close to .42. Part of the discrepancy with 
the model at high latitudes lies in the non-inclusion of a 
surface-ice reflectance parameter. 

To summarize, apart from the region covered by sea-ice, 
there is a marked tendency for the 2/3-CL cloud model to give 
reduced albedos as suggested by Raschke et at. [1973] , Von der 
Haar and Hanson [1969], and others. This result of the 2/3-CL 
model is especially more clearly defined in tropical latitudes 
where Von der Haar and Hanson [1969] have discussed the re- 
ality of both reduced cloud-covers and the resultant reduced 
global -albedo . These latter effects than must give rise to 
larger net radiation into the oceanic surface, which they also 
prescribe . 

In Section VI, heat-budget computations along each of 
the four meridians are presented. These computations make 
use only of the 2/3-CL cloud model for purposes of radiative 
calculations. This selection (2/3-CL) is supported by the 
results of Section VII. However, it is felt that the grid- 
point method used here for 2/3-CL is not the final answer in 
the tropics. In the tropics the peculiar ref lection- charac- 
teristics of cumulus convective elements must be more real- 
istically accounted for than by the large-scale layering 
considered in the solar insolation model of this section. 
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V. SENSIBLE AND LATENT HEAT TRANSPORT AT THE SEA-AIR INTERFACE 



A. GENERAL 

In order to compute heat balances involving the atmosphere 
and the ocean- surface , it was necessary to compute vertical 
fluxes of sensible and of latent heat across the air-sea in- 
terface. The model adopted to describe the turbulent- fluxes 
in the planetary boundary layer was essentially that used in 
the FNWC primitive equation model as delineated by Kesel and 
Winninghoff [1972] and Kaitala [1974] . The model as used in 
this study was adapted in such a way that it gave similar 
quantitative bulk- transfer results across the air-sea inter- 
face as the model used by Kaitala [1974]. This was a major 
deviation from Warner’s study [1974], where he required a 
positive sea-to-air sensible heat transport based primarily 
upon a semi -empirical Bowen ratio. 



B. SENSIBLE HEAT TRANSPORT 

As previously mentioned a slightly different modelling as 
compared to that of Kesel and Winninghoff [1972] was used in 
this study for the vertical distribution of sensible -heat con- 
vergence in the layer 800 to 1000 mb. In the presently opera- 
tional heat-package of FNWC, as described by Kaitala [1974] 
and by Kesel and Winninghoff [1972], the planetary boundary 
layer sensible heat transport is given by 



Hj, (900) 



p 10 C p ( 



K* 



1 + a* ( 9 9 " 6 g 



){ Y, 



e„ - e 



- z 



} 



z 9 " z 10 



10 

(5-1) 
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In (5-1), the parameters are as described below: 



C p = .239 cal gm' 1 (deg. K) ' 1 

0 = T 

x x 

Y c = (the critical countergradient lapse rate) = 

3 x 10 " 5 (deg . K) cm" 1 

K* = (the neutral value of the eddy coefficient) = 

2 x 10 5 cm 2 sec 1 

a* = (a turbulent transfer coefficient) = 

5 x 10 4 cm(deg. K) 1 

P 1Q = (surface density) = p 1Q (R d T 10 ) " 1 • 

The other parameters in (5-1) are 0g , which is the potential 
temperature at 900 mb, 0^ which is the potential temperature 
at the top of the constant flux layer,, and 0 which is the 

O 

potential temperature at the surface. 

Equation 5-1 was given by Kaitala [1974] as the planetary 
boundary layer equivalent of the bulk formula for H^, given 
in terms of surface - layer parameters as 



where is the drag coefficient, V g = surface wind, and 



thin layer between the sea surface and the top of the con- 
stant-flux layer. If Hp of (5-1, 5-2) are always equal, then 
it follows that T x may be solved from the identity, giving 




(5-2) 



(T - T ) is proportional to the temperature lapse for the 
§ X 



T 



x 




(5-3 (a)) 



8b 



where 



K 



1 + 



El 

a* ( 



z 9 




) 



Az = z 9 



10 



R d 9-10 
~g T 




(5-3 (b)) 



(5-3 (c)) 



The requirement of identical values of Hp at level x and at 
level k = 9, with Hp decreasing to zero at k = 8, was inter- 
preted by Kesel and Winninghoff [1972] as giving rise to 
sensible-heat convergence Hp , layer-centered at k = 9 . 
Actually Fig. 10a suggests that the vertical convergence of 
sensible heat flux of amount Hp occurs in the layer (8, 9) 
when the assumption of constant heat flux in the layer (x, 9) 
is made. However, in the sensible-heat model visualized in 
this thesis, a linear decrease in Hp with pressure was as- 
sumed above level 'x', Fig. 10b. Note that level 'x' is as- 
sumed to occur so close to k = 10 that one may solve (5-1) 
for 



Hp(k = 9) = | H 
Solution of (5-4) 



T x = 



C n V T 
D s g 



r (Bulk) . 

by use of (5-1) , 
2 K 0 9 

’ 2 K Y c 

2 K 



Az 



(5-4) 

(5-2) gives the result 
(5-5) 



Equation 5-5 has the same format as (5-3). As was mentioned, 
Kaitala [1974] has suggested the use K* = 2 x 10 5 cm 2 sec 
in (5-3), and this would lead to identical values as the use 
of K* = 10 5 cm 2 sec 1 in (5-5). The latter K* value was tested 
by Warner [1974] who felt, with Kaitala, that either K* value 
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was only an approximate choice of the eddy coefficient re- 
quiring at least some tuning, with a theoretical dependence 
upon the larger-scale synoptic parameters over the ocean. 

The physical implications of condition (5-4) along with 
the linear decrease of Hp(P) with pressure to zero at level 
k = 8 is depicted in Fig. 10b. This condition is compatible 
with uniform convergence of the sensible-heat flux Hp of the 
surface layer into the overlying layer (8, 10). Having uni- 
form convergence from level x (near surface) to level 8 is 
equivalent to maintaining a "well mixed" boundary layer of 
thickness from the surface to 8. 

With either case, Equation 5-3 or 5-5, the sensible heat 
convergences of Fig. 10-a and 10-b were considered applicable 
at level k = 9 and were considered to be identical. In this 
thesis the value of 2 K in (5-5) was based upon an arbitrary 
choice K* = 10 5 cm 2 sec 1 in (5-3b), rather than 2 x 10 5 cm 2 sec 1 . 

In order to evaluate Hp, it was necessary to derive the 
surface winds V g by use of (2-21) at all gridpoints of the 
four oceanic meridians under study. This was done using 
Equations 2-16 and 2-20 with contour gradients extracted from 
the 1000 mb field except in the latitude band ±5°. Here, the 
January climatological values of gradient winds were taken 
from the Tropical Wind Atlas of Atkinson [1970]. 

Besides having computed the applicable surface winds V s 
for use in (5-2), a modification to = 1.4 x 10 (after 
Weiler and Burling, 1967) was introduced, as compared with 
Cp = 2. x 10 3 as suggested by Kaitala [1974]. 
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After determination of T in Equation 5-5, Hp was then 
computed at each gridpoint using Equation 5-2. The values 
of Hp are used later in calculations of heat budgets in 
Sections VI and VIII. 

C. EVAPORATIVE HEAT TRANSPORT 

As in the case of sensible heat convergence, the latent 
heat transport by turbulence was assumed to be subject to 
uniform heat flux convergence throughout the layer 800-1000 mb 
(Fig. lib). The total amount of latent heat removed by evapo- 
ration from the ocean surface was similar to that given by 
Kaitala [1974] using bulk transfer theory. While a turbulent 
model similar to Fig. lib can be adopted, it is necessary to 
recognize that the realized latent heat may occur at arbitrary 
levels (for example, higher than k = 8) and therefore is not 
introduced at level k = 9. The basis of the turbulent latent- 
heat flux model adopted here is shown in Fig. lib, where 



E (goo) = %E (Bulk) (5-6) 

and where the surface - layer bulk transfer of latent heat is 
given by 



E = E(Bulk)=Lp 10 C D V s (q 1Q -q x ). 



(5-7) 



From the planetary boundary- layer transfer theory analog to 
Hp(900), E(900) may be written 



E (900) Lp 10 K z 



~ ^ q x' q 9^ 



(5-8) 



'9 "10 

In (5-7), and (5-8), L is the latent heat of vaporization 



L = 596.73 - .601TpQ cal gm , T 1n in °C. 



10 
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Also in (5-8) and q x are mixing ratios at the level k = 9 

and at x, the top of the constant flux layer. The factor h 

which appears in the right side of (5-6) gives the type of 

latent heat convergence as shown in Fig. lib, as against 

that of Kaitala (Fig. 11a). Combining Equations 5-6, 5-7, 

and 5-8 leads to the solution for q , namely 

X 



*x = 



C D V s q 10 + 2 Az q < 
C D V s + 2 h 



(5-9) 



with Az given by (5-3(c)). As in the case of T in Equation 
(5-5), the factor of 2 K of (5-9) was taken equal to Kaitala's 
K of Equation 5-3(b). This was done in this study by using 
a value of K* = 1. x 10 5 cm 2 sec 1 as against 2 x 10 5 cm 2 sec 1 
suggested by Kaitala. Thus amounts computed for E of Equa- 
tion (5-7) were equivalent to amounts computed using Kaitala's 
scheme (Equation 5-8) only the turbulent convergence of E 
was distributed differently in the planetary boundary layer. 
Thus after computing q (Equation 5-9), E was calculated at 
each gridpoint except those over ice by using Equation (5-7). 



D. TURBULENT HEAT TRANSPORTS OVER AN ICE -COVERED OCEAN 

Any gridpoint whose surface temperature was less than -5°C 
was considered to be over an ice-covered ocean for the purposes 
of calculating H^ and E. In this study there were four such 
gridpoints. In the handling of H^, a formula adopted from 
Kaitala [1974] was first utilized to describe heat conduction 
into the ice (HICE) , from the warmer water beneath the ice 
surface : 



HICE = BBB* (T - T ) 

o g J 



(5-10) 



where BBB is an ice conduction coefficient (697.83 ergs 

- 2 “1 - l 

cm sec (deg. K) ) , T q is the freezing point of sea water 
(271.3 deg. K) , and T is the surface temperature. HICE is 

o 

then utilized in the following surface heat-balance equation 
to determine the summation of the sensible heat flux Hp plus 
the latent heat flux: 

H r + E = HICE + QABG - F Q *. (5-11) 

A Bowen ratio r as described by Kaitala [1974] was then cal- 
culated using Equation 5-12, where the right side is an em- 
pirical least-squares best-fit applicable over continental 
and ice-covered ocean areas. 

Hp 

r(4>) = — g- = 9 . 6 | Sin 4> | 2 - 7.93|Sin <f> | + 2.0. (5-12) 

Here cf> is the latitude of the gridpoint, and r was never al- 
lowed to exceed 1.33. This restraint on r meant that E was 
always bounded by 

E < .75 Hp (5-13) 

(cf . Kaitala, 1974) . From the definition of the Bowen ratio 
(r) as a function of Sin <j>, and using the following equation 
for the combined interface fluxes Hp + E as 

Hp + E = Hp (1 + ~ ). (5-14) 

one obtains simply 

H r = ttt ( H r + E) * (5-iS) 

Since both the Bowen ratio (r) and Hp + E are known at each 
gridpoint, Hp is easily determined from (5-14). Then a simple 
rearrangement of terms in the Bowen ratio definition permitted 
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the computation of E. This is given by 



E 



H, 



r (40 



(5-16) 



Equations 5-10 and 5-11 always gave negative or downward 
values of sensible and evaporative heat fluxes at the four 
gridpoints involving ice-covered oceans since QABG - F^o* 
was dominantly negative. Further, evaporation never exceeded 
75% of Hp in magnitude. These negative values appear to agree 
generally with those arrived at in the work of Vowinckel and 
Taylor [1964] which dealt in part with evaporative and sen- 
sible heat fluxes over an ice-covered Polar Ocean areas in 
winter . 



E. LARGE-SCALE TURBULENT HEAT FLUX ACROSS THE AIR-SEA INTERFACE 
At each gridpoint, it is possible through the adaptation 
of the turbulent boundary layer model to include heating 
rates across the air-sea interface. For use in Section VI, 
we compute the values of Hp as the convergence of sensible- 
heating rate for the layer (8, 10) through Equations 5-2,..., 
5-5. Also for Section VI, the combination E + Hp considered 
primarily as a heat loss at the ocean surface is computed at 
each gridpoint by use of Equations 5-7, 5-9. For the over- 
ice cases Hp + E is determined by (5-11), while Hp results 
from (5-12) and (5-14). In any case, the model governing 
these heating rates, Hp and E has been extracted from the FNWC 
operational turbulent -boundary layer model and appended to the 
radiational model (with the 2/3-CL case) for operational com- 
putations at each gridpoint of the four meridians of interest. 



93 



The values of Hp turned out to be in general small and 

negative in low latitudes reflecting' computed negative values 

of T - T . However, in the data-rich mid-latitude cross 

g x 

sections where both strong surface winds and relatively warm 
ocean temperatures coexist, Hp is positive and larger and 
dominates the net atmospheric heating in layer (8, 10). Re- 
gardless of the inconsistency in the sign of Hp , the values 
of E + Hp seems to be reasonable heat-loss rates at the ocean 
surface at all latitudes insofar as comparison with other 
sources of data (Budyko, 1956) are relevant. 

Finally in the large-scale sense Hp + E is included as a 
(positive) contribution to the heat budget of the tropospheric 
column even though it is not known precisely at what level the 
latent heat flux E will be realized. 
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VI. MERIDIONAL CROSS-SECTIONAL DEPICTION OF THE HEAT BUDGET 
SYSTEM COMPUTATIONS 



A. GENERAL 

The general design of this section is to utilize all com- 
putational subroutines explained in Sections III, IV and V 
in computations of a single time-step of the heating model 
developed for use in the FNWC prediction model. It was de- 
cided after the testing of F£ and of ALB with the correspond- 
ing values of Raschke et al . [1973] that only the two-thirds 

cloud model would be displayed in the meridional cross- 
sections depicting the single time-step computations of the 
heating package of the FNWC prediction system (Figs. 12, 

13(a), 13(b), 14(a), 14(b), 15(a), 15(b), 16(a), and 16(b)). 
The appropriate radiative calculations with two-thirds CL 
were thus performed at each gridpoint of the four meridians 
used for presentation purposes in this section. However, the 
same type computations were made with full-cloud cover so 
that the mean radiative balances could be compared layer-by- 
layer and at the levels k = 2 and k = 10. Both sets of com- 
putations were zonally-averaged and used to get mean radia- 
tion meridional cross-sections displayed in Section VII 
(Figs. 17, 18(a), 18(b), 19(a), and 19(b)). 

In this study, the FNWC gridpoint processed analyses for 
0000GMT, 16 January 1974 were used at the three Pacific cross- 
sections, while that for 1200GMT, 16 January 1974, was used 
for the single Atlantic meridian. In any case, the set of 
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turbulent transfer values are in ly min 1 for the levels or layers considered. The 
reduced cloud amounts 2/3 CL(1) and 2/3 CL(2) are given by Eq. (2-15). 



meridians was considered as representative of 16 January 
data when converted to heating-rate calculations. 

Figure 12 depicts in symbolic language the key to the 
computational entries in Figs. 13, 14, 15, and 16. In 
actuality, this key list represents computations made at 
each (I, J) gridpoint sounding from data in the form of 
Table 1(b). It covers those computations proceeding from 
the top of the troposphere to the ocean surface, including 
sensible and latent heat transfers (Equations 5-2, 5-7). 

For simplicity and ease of climatological data comparison, 
Figs. 13, 14, 15, and 16 were developed by interpolating 
gridpoint results to whole 5° latitude increments. The in- 
terpolation routine utilized was the Lagrangian cubic inter- 
polation scheme. It has the form 



where I is determined from an equation similar to (4-3(c)), 

Q(I) is the desired interpolated value of the heating parameter 
at a non-integral grid-value I that corresponds to a whole 
latitude of 5°. In (6-1), the indices 1, 2, 3, and 4 are 
for any four increasing successive integral values of I along 
gridpoint lines. The symbol I in (6-1) is usually applied 
between points 2 and 3, but this interpolation scheme may 
also be applied with good accuracy to an I-value which is 
non-centered in the I-array (1, 2, 3, 4). In order to com- 
plete the interpolation to latitude, the array is applied 
sequentially to successive sets of integral I-values along each 



Q(I) = Q 



(1-2) (1-3) (1-4) 
1 (1-2) (1-3) (1-4) 




+ o (I-D(I-2)(I-4) + (1-1) (1-2) (1-3) 

^3 (3-1) (3-2) (3-4) (4-1) (4-2) (4-3) 



meridian. 



97 



B. EXPLANATION OF TERMS IN FIGURE 12 



1 . Cross-Section at Level k = 2 

It should be explained here that all insolational 
parameters enumerated upon in Section IV dealt with the speci- 
fic time of day that corresponded to the hour angle h at the 
instantaneous times under consideration. The resultant in- 
cident solar insolation dealt with is then by Equation 4-1 

F ( 2) = S (“— ) Cos z. 
m 

In order to avoid reference to specific map times, the applic- 
able solar hour-angles were: h = 35°, 10°, 55°, and 35°, 

respectively, for cross-sections 1, 2, 3, 4, depicted in 
Figs. 13, 14, 15, and 16. 

Regarding the first symbol at level k = 2 in Fig. 12, 
note QAVE and not F(2) is displayed. This is because averages 
over the 24-hour day (16 January 1974) are considered more 
meaningful for the data day under consideration, and hence 
are presented in this form in Figs. 13, 14, 15, and 16. QAVE 
is the 24 hour average of F ( 2 ) found by using a similar form 
of Equation 4-1 but replacing Cos z by Cos z as in (6-2) 

QAVE = F (2) ( §|§-| ) (6-2) 

where 

Cos z = [H Sin 4> Sin 6 + Cos 4> Cos 6 Sin H] (6-3) 

IT 

H(hour angle at sunset) = arc Cos[- Tan 4) Tan 6] 

<f> = latitude 

6 = solar declination for 16 January 1974. 
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Cos z is equal to the 24-hour average cosine of the zenith 
angle. The 24-hour time-averaging period is considered con- 
sistent with the longwave flux divergences, which change 
slowly with the time of day. This conversion to expected 
daily averaged heat-transfer quantities is compatible with 
the determination of a heat budget for the given date (e.g., 
16 January 1974) . 

In keeping with the concept of using 24-hour averages 
other parameters for level k = 2 are defined as follows: 

QREF = REF(t) ( ) (6-4) 

where 

REF (t) = F (2) - A26 - A610 - (IA10 + IS10) . (6-5) 

Here REF(t) is the instantaneous solar- ref lected insolation 
at the gridpoint as given in 4-30 and QREF is its 24-hour 
average, assuming that planetary albedo remains constant for 
the 24-hour period, even though it was calculated for the in- 
dicated solar time at each meridian. 

The same principle will be used with regard to all 
other solar disposition parameters in converting from time- 
dependent values at solar time t to 24-hour averaged values. 
Thus in the following Equations 6-7, 6-9 below Q26(t) be- 
comes Q26, and Q610(t) becomes Q610. Superior bars ( ) will 
not be explicitly shown in the values presented in the cross- 
sections key, Fig. 12, but will be understood in what follows 
Finally, the 24-hour average "system" balance BALT is com- 
puted from 

BALT = QAVE - (QREF + F 2 *) (6-6) 
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for the tropopause level k = 2 at the indicated latitude. In 
computing this balance the value of F£* was introduced from 
Equation 3-11. Net terrestrial fluxes "F*" such as F£* that 
appear in Fig. 12 were considered to be constant throughout 
the 24-hour period which would be approximately true if the 
cloud-covers were representative of the period. 

2 . Cross-Section in Layer (2, 6) 

Referring to Fig. 12 the following definitions apply. 

In this layer all the heat-transfers shown are assumed to 
be radiative only in our version of the FNWC heating model. 

The averaged radiative cooling (heating) rate is given by 

BAL26 = Q26 - F26 (6-7) 

where 

Q26 r daily average absorption in layer (2,6) and is 
defined relative to A26(t) by a transformation 
similar to Equation 6-2 

F26 = terrestrial cooling rate given by Equation 
3-12 (a). 

3 . Cross-Section in Layer (6, 10) 

When considering radiative transfers only in the layer 
(6, 10) the 24-hour average radiative cooling is given by 

BAL610 = Q610 - F610 (6-8) 

where Q610 is the sum of Q68 and Q810 in Fig. 12. BAL68 has 
been taken as one-half of BAL610 . However 

BAL810 = Q810 - F810 + H f = BAL68 + Hj, . (6-9) 

Here Hp is the sensible heat flux transported to (8, 10) by 
Equation 5-2. Likewise, F610 is the sum of F68 and F810, 
terrestrial flux loss values for the respective layers. Note 
here that the parameter BAL610 remained negative at all gridpoints 
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of Figures 13 through 16 when considering radiative transfers 

only. 

4 . Cross-Section at Air-Sea Interface (k = 10) 

The heat balance at the earth's surface (BALB) con- 
sisted of the following equation 

BALB = QABG - F 1() * - (Hp + E) . (6-10) 

Here QABG is the 24-hour average insolation absorbed by the 
surface as 

QABG = QABG(t) [Cos z / Cos z] (6-11) 

(E + Hp) is the sum of the turbulent heat transfers 
which were computed from the bulk formulas as explained in 
Section V.B. and V.C. 

C. MERIDIONAL CROSS-SECTIONS OF THE VERTICAL HEAT BUDGET 
FOR 16 JANUARY 1974 

Figures 13, 14, 15, and 16 as explained represent the 
single-time step of heating computations for each of the four 
meridians used in this study. Note that the four Figures 
(13-16) are divided into a and b parts with section (a) 
showing "tropical" results and (b) mid-to-high-latitude re- 
sults. While the results depicted in these cross-sections 
are exhibited as representing daily - averaged values, they are 
actually based upon heat-budget parameters at the specific 
map time 0000GMT and/or 1200GMT on 16 January 1974. Hence, 
if the results are to be used as part of the "heat-package" 
subroutine of FNWC, all solar-radiative absorption and-reflec- 
tance terms must be recovered as functions of GMT e.g., 

F ( 2 , t) = QAVE * [Cos z / (Cos z)] 
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and 



QREF(t) = QREF * [Cos z /(Cos z) ] 

etc. These solar disposition terms may then be utilized in 
connection with a one-hour time-step application of the 
thermodynamic equation of the set of primitive equations used 
in the FNWC prediction-process. The purpose for leaving all 
results of Figures 13,..., 16 in the form of daily-averages 
was primarily so that these figures could be compared with 
climatology where available of the zonally-averaged heat 
budget of the ocean-atmosphere system for January. 

Finally while only the (2/3) -CL cloud model results have 
been depicted in Figs. 13, 14, 15,16. An analogous set of 
heat-budget cross-sections were developed for the full-CL 
model. These are not presented here because of the sheer 
volume of the results. However, certain radiative terms in 
both the (2/3) -CL and full-CL cross-section sets have been 
extracted from both sets for purposes of presenting zonal 
comparisons of the purely radiative-heat -budget of the FNWC 
data as it was applicable to 16 January 1974. The compara- 
tive zonally-averaged radiative cross-sections are presented 
in Figs. 18 and 19 of Section VII. 
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VII. ZONALLY-AVERAGED RADIATIONAL BALANCE OF THE OCEAN- 
ATMOSPHERE SYSTEM 



A. GENERAL 

In order to depict the results of the purely radiative 
contributions, zonally-averaged over the ocean-atmosphere 
system for this study, Figs. 18(a, b) and 19(a, b) are pre- 
sented. These figures present respectively the zonally- 
averaged radiative cross-sections for both the 2/3-CL and 
the full-CL cloud amounts. These cross-sections show the 
results after zonal -averaging over the four meridians con- 
sidered in this study. The results are then displayed in 
the format of Fig. 17, which gives the key to the line items 
of the cross-sections. Note that Fig. 17 simply combines 
the radiative- transfers Q68, Q810, and also F68 with F810 
of Fig. 12 to arrive at layer- average radiational warming 
rates in the combined layer (6, 10). Also in the results 
of this section, turbulent - trans fer heating rates into the 
atmosphere (and corresponding cooling in the underlying ocean) 
have merely been omitted from the results of Figs. 13,..., 16, 
in passing to the radiative composites presented here. 

In obtaining a simple zonally-averaged set of radiative 
heating-rates at each latitude <f> , all values of each parameter 
listed in Fig. 17 at latitude <J) in the range 20S,...,65N (by 
5° steps) were simply averaged over the meridians providing 
the averaged value at <f> . It should be noted that at <j> = 65N, 
there was only one contribution to the average. In the Southern 
Hemisphere region 20S,...,10S, only two values (from 125W and 
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LAT latitude of the gridpoints in 5° increments 

a) QAVE 24-hour averaged mean insolation at level k=2 

b) QREF mean reflected average insolation at level k=2 

c) F 2 * mean net outgoing long-wave flux at level k=2 

d) BALT mean averaged earth tropospheric gain rate (a-b-c) 
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Figure 18(a). Zonally-averaged radiational cross-section for tropical latitudes with 
2/3-CL model. Refer to Fig. 17 for key. All values listed are daily averages in 
ly min -1 and are computed from data for 16 January 1974. 
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Figure 18(b). Zonally-averaged radiational cross-section for mid to high latitudes 
with 2/3 CL model. Refer to Fig. 17 for key. All values listed are daily averages 
in ly min -1 and are from data for 16 January 1974. 
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all radiative calculations. 
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Ea. 2-14) is used in all radiative calculations. 



and 35W) of each parameter contributed to the means listed 
here in Figs. 18 and 19. At all other latitudes, there were 
in general four values of each radiative parameter entering 
into the computed radiative mean cross-section. As a re- 
sult, one may conclude that near the northern and southern 
boundaries of Figs. 18 and 19, the listed values of the 
zonally-averaged quantities may not be exact, but the general 
zonal trend is probably reliable. (Fig. 20 shows this to be 
a valid result.) 

Figures 18, 19, and 20 are useful for the discussion of 
the earth-atmosphere system radiative balance , and in the 
discussion of the relative merits of the two-thirds cloud- 
model as compared with the full-CL model. The full-cloud 
model used by Warner [1974], is also used in comparative 
discussion of the heat budget of the earth-atmosphere system 
of Section VIII. 

B. EARTH -ATMOSPHERE SYSTEM RADIATIONAL BALANCE RESULTS 

Figure 20 depicts the zonally-averaged distribution of 
three parameters, R , R , and R for full cloud amounts and 
two-thirds cloud amounts, respectively. These parameters 
were defined (after Malkus , 1962) as (1) R g , the mean radia- 
tive energy transfer rate across the top of the ocean-tropo- 
spheric system (referred to as BAL'T at k = 2 in Fig. 17; (2) 
R o , the mean radiative cooling rate in the troposphere, i.e., 
the sum of BAL26 and BAL610 in Fig. 17; and (3) R, the mean 
radiative warming (cooling) rate at the earth's surface, 
referred to as BAL10 in Fig. 17. The relationship bet\\reen 
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Figure 20. Radiational balance at the tropopause (R ), at 

s 

the ocean surface (R), and in the tropospheric column (R ). 

cl 

Solid lines denote computations made with the 2/3-CL (sub- 
scripted "1") cloud model and dashed lines correspond to the 
full-CL model, (subscripted "2"). All computations are for 
16 January 1974 on ly(min) x 10 2 . 
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these three radiative parameters is as follows: 

R s = R + R a . (7-1) 

Table V lists zonally-averaged comparisons of R , R , and R 
as a function of latitude <j> for this study for 16 January 
1974, both for the 2/3-CL and the full-CL cases. Also in- 
cluded is the R g ^ values (signifying satellite climatology) 
after Raschke et al. [1973], and R 7 and R x values. The 
latter two symbols denote December climatological values 
after Berliand [1956] corresponding to R and R respectively. 
"Hemispheric" means of any Q-value are then by the following 
cosine -weighting scheme: 

65 

• E n Q. Cos <J>. 

0 = — ^ i i- f7-21 

V Wt.Avg. 65 * ^ 

. Cos <}> . 

1 = 0 r l 

The tabular values of R , R, and R as a function of 

s’’ a 

latitude and for the two cases, 2/3-CL, and full-CL, have 
been plotted in Fig. 20, where some interesting results of 
the radiational model become evident. 

It is evident from Fig. 20 that the R values obtained 
with the 2/3-CL cloud model exceed those based upon the full- 
CL model up to 4> = 40N. Poleward of this latitude the two 
R g -curve cross with R g (CL) becoming slightly larger (although 
negative) than R (2/3-CL) . The same general results hold for 
both R^ and R£ and for R ^ relative to R up to latitudes 
40-45N, where also cross-overs in the graphical depictions 
occur in both cases. These results are related to the 
zonally-averaged radiative budgetary cross-sections of Figs. 
18 and 19 as is explained in Section VII (C) . 
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Earth-Troposphere 
Net Radiation, R 
(k=2) s 


Surface Net 
Net Radiation 
(k=10) 


, R 


Atmospheric-Net 
Radiation Gain 
R 

a 


Lat 


R 

S 4 


R 

s i 


R 

S 2 


Ri 


r 2 


r 3 


R 

a i 


R 

a 2 


R 

3 3 


20S 


.1800 


.0681 


.0216 


.1622 


.1123 




-.0942 


-.0913 




15 


.1700 


.0552 


.0052 


. 1624 


.1125 




-.1073 


-.1075 




10 


.1600 


.0488 


-.0039 


.1780 


.1310 




-.1292 


-.1347 




5 


.1430 


.0475 


.0046 


.1902 


.1537 




-.1426 


-.1492 




0 


.1150 


.0182 


-.0270 


.1802 


.1443 


.1713 


-.1621 


-.1720 


-.0833 


5 


.0825 


-.0050 


-.0478 


.1697 


.1387 




-.1746 


-.1864 




10 


.0550 


-.0140 


-.0428 


.1645 


.1454 


.1481 


• -.1787 


-.1882 


-.1157 


15 


.0150 


-.0489 


-.0736 


.1374 


.1226 




-.1863 


-.1962 




20 


-.0125 


-.0813 


-.1010 


.1123 


.1022 


.1088 


-.1937 


-.2031 


-.1551 


25 


-.0256 


-.0938 


-.1058 


.0807 


.0738 




-.1745 


-.1796 




30 


-.0400 


-.0998 


-.1041 


.0475 


.0429 


.0463 


-.1472 


-.1471 


-.1551 


35 


-.1050 


-.1497 


-.1551 


.0152 


.0145 




-.1648 


-.1696 




40 


-.1375 


-.1689 


-.1706 


-.0092 


-.0072 - 


.0116 


-.1597 


-.1634 


-.1597 


45 


-.1625 


-.1884 


-.1835 


-.0094 


-.0068 




-.1791 


-.1768 




50 


-.1750 


-.2062 


-.1969 


-.0275 


-.0200 - 


.0417 


-.1788 


-.1769 


-.1828 


55 


-.2025 


-.1798 


-.1671 


-.0409 


-.0300 




-.1388 


-.1371 




60 


-.2300 


-.1932 


-.1759 


-.0868 


-.0683 - 


.0509 


-.1065 


-.1077 


-.2037 


65N 


-.2400 


-.2265 


-.2265 


-.0868 


-.0868 




-.1397 


-.1397 




Wt. 

Avg 

(N. 

Hera) 


-.0500 


-.1000 


-.1139 


.0668 


.0581 


.0699 


-.1668 


-.1720 


-.1437 



Table V. Zonally- averaged components R, R , R of the earth-atmosphere 

s a 

system as a function of latitude. The averaged values as found in this 

study using subscript "l" for the 2/3-CL case and the subscript ,f 2" for 

the full-CL case. Corresponding values of R^ reported by Raschke et al. 

(1973) are denoted by subscript "4 M and those of R and R reported after 

a 

Berliand (1956) for December climatology denoted by subscript "3". 
Weighted averages over the latitude range 0-65N are shown in the last 
line of the table. (All values in ly min” 1 ). 
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C. EXPLANATION OF THE MEAN ZONAL RADIATIVE - BUDGET DISTRIBUTIONS 



It should be observed that in Table V, the R values ob- 
tained with the 2/3-CL model are greater than those resulting 
from the full-CL model throughout the tropics and up to lati- 
tude 40N (approximately) . Poleward of latitude 40N, the 
values of R s (2/3-CL) become smaller than those for full-CL. 
These comparisons are illustrated clearly in Fig. 20 by the 
zonal depictions of the R g values calculated by the two 
cloud-model cases. Examination of Table V for comparative 
values of the radiative balances at the surface shows the 
same general effect at the ocean surface, with the cross-over 
between 2/3-CL and full-CL cases (R-^ and R£ curves in Fig. 20) 
occurring near latitude 40N. Finally, the same general effect 
is evident in the comparative values of R a (2/3-CL) compared 
to R (CL) . That these three comparisons should follow each 
other systematically is perhaps to be expected in view of 
the connecting relationship 

R = R + R . 
s a 

Of the changes resulting from substituting (2/3-CL) for the 
full-CL cases, the largest gain in low and mid- latitudes oc- 
curs in the comparative R g -curve for the (2/3-CL) case. The 
smallest gain occurs for the computed R^-changes, and an in- 
termediate gain is registered for the R(2/3-CL), the net 
radiative flux at the ocean surface. 

The same results for R s -changes may be seen by comparing 
Fig. 18 for the 2/3-CL case with Fig. 19 for the full-CL case. 
In such a comparison it is seen that the systematic reduction 
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of CL to 2/3-CL increases the net useful solar insolation 
at k = 2, (Figs. 17, 18, and 19) by reducing QREF. However, 
the reduction achieved in QREF(2/3-CL) is quite small at 
latitudes 4> > 45N. On the other hand with the 2/3-CL case, 
the F 2 * values (net terrestrial loss to space at k = 2) is 
increased by a small but consistent amount which at higher 
latitudes <f> > 45N, exceeds the very small net solar insola- 
tion gain (k = 2) . 

It therefore turns out that in the consideration of 
Figs. 18 and 19 of the atmospheric- layer radiative transac- 
tions (both solar and terrestrial) that the net radiation 
R penetrating to the surface is influenced primarily by QREF, 
at best in lower latitudes and up to <j> = 40N. Here -A (QREF) 
exceeds A as a result of the change in cloud cover and 

thus accounts for greater R g -values. The same effects ac- 
count for R-changes at the surface, while at <t> > 45N the 
cross-over effect holds both at k = 2, and also at k = 10. 

The curves of R versus 4> in Fig. 20, in effect measure 
the atmospheric participation in the radiative absorption and 
emission processes, follow the lead of -A(QREF) in lower 
latitudes and of F£* in higher latitudes, which act as the 
forcing functions in the radiative balance of the system. 

Recent satellite results especially of Von der Haar and 
Hanson [1969] as well as of Raschke et at. [1973], also draw 
attention to the reduced values of global albedo, compared 
to previous heat budget estimates (e.g., London, 1957) in 
tropical and subtropical latitudes. This also shows up in 
Table V in the comparative values of Raschke 's distribution 
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of R g with the "best" R s distribution computed here for the 

2/3-CL model. It should be noted from Table V that the use 

of the 2/3-CL model leads to R , -values which are in closer 

s 1 

agreement with R^ up to latitude 40N, than that generated 
by the full-CL model. At high latitudes in Table V the 
comparisons between the results of these two cloud models 
are inconclusive in their relationship to Raschke's results. 
However, his results in some of these areas are over ice, 
which enhances his effective QREF as compared with either 
model used here. 

Finally note that comparative values after Berliand 
[1956] are also included in Table V for R and R at 10° lati- 
tude intervals based on December global climatology. The 
most noteworthy comparisons to be made here is the transi- 
tion from positive to negative surface net flux R- which oc- 
curs between 35-40N, in reasonable agreement with the results 
interpolated either from the R-^-, R 2 "distribution of Table V. 
In general, the global results of Berliand reflect a greater 
continental trend with latitude than is found in the 
January 1974 study over the oceans. 



123 



VIII. ZONALLY -AVERAGED TROPOSPHERIC AND OCEANIC HEAT BUDGETS 
FOR 16 JANUARY 1974 



A. THE TROPOSPHERIC HEAT BUDGET OVER ICE -FREE OCEANS 

By averaging over the tropospheric column for the four 
meridians displayed in Figs. 13,..., 16, latitude by latitude, 
the zonally-averaged tropospheric balance results (Table VI 
and Fig. 21). Values presented in Table VI and Fig. 21 are 
R , the tropospheric radiative net cooling defined in Section 
VIII. B., and E + Hp as formulated in Equations 5-2, 5-7. The 
latter parameter represents the combined turbulent heat-trans- 
fer lip + E across the ocean-air interface (see for example, 
Figs. 13,..., 16). The heat-budget equation for the atmosphere 
column at latitude <f> then becomes, with the notation of 
Malkus [1962]: 

Q V a * S a * R a + + H r> • f 8 ' 1 ) 

Here S is the in-situ storage-heating rate of the column, 

3 

and Q va is the required flux-divergence of heat necessary to 

bring about the "balance" in the column, for the observed 

values of R E + H„ and S . 

a I a 

Table VI lists the zonally-averaged values of R , E + Hp 
and the resulting values of Q va + S , which are computed from 
(8-1). As before the subscript notation "1" and "2" denote 
the 2/3-CL and full-CL computations respectively. Figure 
21 depicts each function R , E + Hp and Q va + S a , computed 
by the two cloud-models. 

As noted before in Section VII. C., R tends to exceed 

’ al 

R Q 2 over most of the tropical Southern Hemisphere and in the 
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Tropospheric 
Radiative heating 
rate R 



Turbulent transport 
heating rate in 
troposphere 



Tropospheric Balance 

requirements 

Q +S 
va a 



Lat 



(E+H ) (E+H p ) (E+H r ) 
1 1 1 2 i 



(Q +S ) (Q +S ) (Q +S ) 
va a i va a 2 va a 



20S 


-.0942 


-.0913 




.1402 


15 


-.1073 


-.1075 




.1718 


10 


-.1292 


-.1347 




.1689 


5 


-.1426 


0.1492 




.1197 


0 


-.1621 


-.1720 


-.0833 


.0839 


5 


-.1746 


-.1864 




.0791 


10 


-.1787 


-.1882 


-.1157 


.1374 


15 


-.1863 


-.1962 




.1611 


20 


-.1937 


-.2031 


-.1551 


.1430 


25 


-.1745 


-.1796 




.1257 


30 


-.1472 


-.1471 


-.1551 


.1998 


35 


-.1648 


-.1696 




.2341 


40 


-.1597 


-.1634 


-.1597 


.2971 


45 


-.1791 


-.1768 




.2172 


50 


-.1788 


-.1769 


-.1828 


.1901 


55 


-.1388 


-.1371 




.1301 


60 


-.1065 


-.1077 


-.2037 


.1156 


65N 


-.1397 


-.1397 




-.0721 


Wt. 

Avg 

(N. 


-.1668 


-.1720 


-.1437 


.1508 



Hem) 



.1402 




.0461 


.0490 




.1718 




.0644 


.0642 




.1689 




.0397 


.0342 




.1197 




-.0229 


-.0294 




.0839 


.0764 


-.0782 


-.0882 


-.0069 


.0791 




-.0956 


-.1073 




.1374 


.1018 


-.0413 


-.0508 


-.0139 


.1611 




-.0252 


-.0351 




.1430 


.1342 


-.0507 


-.0601 


-.0209 


.1257 




-.0488 


-.0538 




.1998 


.1250 


.0526 


.0528 


-.0301 


.2341 




.0692 


. 0644 




.2971 


.1134 


.1374 


.1337 


-.0463 


.2172 




.0382 


.0404 




.1901 


.1458 


.0113 


.0132 


-.0370 


.1301 




-.0087 


-.0070 




.1253 


.1944 


.0091 


.0172 


-.0093 


.0721 




-.2118 


-.2118 




.1512 


.1208 


-.0160 


-.0208 


-.0229 



3 



Table VI. Zonally-averaged components R , (E+H p ) and (Q +S ) of the 

a I va a 

tropospheric heat budget. The averaged values as found in this study 
using subscript "l" for the 2/3-CL case and the subscript "2" for the 
full-CL case. Corresponding values reported after Berliand (1956) for 
December climatology denoted by subscript u 3". Weighted averages over 
the latitude range 0-65N are shown in the last line of the table. (All 
values are in ly min” 1 ). 
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Figure 21. Distribution of R„ , E+H r and of Q, +S o as a 

cl 1 V d ci 

function of latitude for the two cloud-model cases. 

(solid curves — 2/3-CL case; dashed curves -- full-CL case) 
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Northern Hemisphere to close to 40N, although these differences 

are quite small. The distribution of E + as a function of 

latitude is independent of the cloud-model except at high 

latitudes over a frozen sea-surface which occurs near <f> > 60N. 

Figure 21 shows the identical sets of E + up to latitudes 

55N. The curves of Q + S are also shown in Fig. 21, as 

va a 

the residual of the right sides of Equation 8-1 for both the 

2/3-CL and full-CL cases, respectively. Note that (Q va + S a ) ^ 

slightly exceeds (Q + S )~ over most of the latitude range 

v a a l 

to 25N as a result of the corresponding of excess R . over 

a -L 

R & 2 for the same latitude range. Note also that both 
Q va + S a distributions reach peak values at <p = 40N, and 
that the values Q va + S a are generally positive over 30-50N. 

This latter feature follows from (8-1), since Q + S 
must follow the trend of E + where the latitudinal varia- 
tion of the latter function peaks sharply, relative to more 
or less uniform behavior of R (4>) . 

The hemispheric mean value of XO + S ) = -.0160 ly min 1 

V ci 3 

per cm 2 column. In this study, Q, r „ has not been computed, but 
if the XO + S ) is attributed to mean storage-cooling S 
of the troposphere alone, the daily storage rate corresponds 
to a temperature -change rate given by 

w (R~ + E + HjO 

( It V 4 - 1 — JFH E — - x 1440 min/da ^ (8 ' 2) 

with AP mb = 800 mb in the troposphere. This resultant mean 
in-situ hemisphere cooling rate (8-2) over the tropospheric 
depth of 800 mb is then 
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it ) a = - . 12°C (day ) 

averaged over the cm 2 tropospheric column. This seems to be 
a reasonable midwinter global value. 

The variation in the computation of E + in the 2/3-CL 
and the full-CL cases which occurs poleward of latitude 55N 
is explained in Section VIII. C. 

B. THE ZONALLY-AVERAGED HEAT BUDGET OF THE ICE- FREE OCEAN 
For the present, neglect the minor variations attributed 
to water-to-ice conduction rates which occur in latitudes 
poleward of 55N (see for example H. of Equation 5-10). 
Neglecting this ice -conduction term, the zonally-averaged 
heat-budget of the ocean water-mass obeys the thermodynamic 
relationship 

Q vo + S = R + [ - (E + H p )] (8-3) 

for both the 2/3-CL and full-CL cases. Note that Equation 
8-3 is the analog of (8-1) for the atmosphere, with -(E+H^) 
representing generally the cooling rate for the water mass. 

The terms of the right side of (8-3) have been computed 
for each meridian and averaged by latitude, with the results 
shown in Table VII. Then Q vq + S is computed as the residual 
of the sum of the two functions on the right side of (8-3), 
and shown also in columns 1 and 2 of Table VII. The graphi- 
cal distribution of the results listed in this table are 
shown in Fig. 22 both for the 2/3-CL (subscript "1") and 
the full-CL (subscript "2") cases. 

As before the primary ocean-mass heating- function is R, 
the radiative heating rate at the ocean-air interface. For 
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Surface Net 
Radiation, R 
(k=10) 




Turbulent transport 
heating rate loss 


Surface Balance 
Requirements 


Lat 


Ri 


r 2 


Ra 


“(E+Hj,) i 


- (E+H p ) 2 


- (E+H p ) 3 


(Q vo +S) 


1 (Q vo +s) 2 


^vo +S) 3 


20S 


.1622 


.1123 




-.1402 


-.1402 




.0220 


-.0280 


15 


.1624 


.1125 




-.1718 


-.1718 




-.0094 


-.0592 




10 


.1780 


.1310 




-.1689 


-.1689 




.0091 


-.0379 




5 


.1902 


.1537 




-.1197 


-.1197 




.0705 


.0340 




0 


.1802 


.1443 


.1713 


-.0839 


-.0839 


-.0764 


.0964 


.0604 


.0949 


5 


.1697 


.1387 




-.0791 


-.0791 




.0906 


.0596 




10 


.1645 


.1454 


.1481 


-.1374 


-.1374 


-.1018 


.0271 


.0080 


.0463 


15 


.1374 


.1226 




-.1611 


-.1611 




-.0238 


-.0385 




20 


.1123 


.1022 


.1088 


-.1430 


-.1430 


-.1342 


-.0306 


-.0408 


-.0254 


25 


.0807 


.0738 




-.1257 


-.1257 




-.0450 


-.0520 




30 


.0475 


.0429 


.0463 


-.1998 


-.1998 


-.1250 


-.1524 


-.1570 


-.0787 


35 


.0152 


.0145 




-.2341 


-.2341 




-.2188 


-.2195 




40 


-.0092 


-.0072 - 


.0116 


-.2971 


-.2971 


-.1134 


-.3063 


-.3043 


-.1250 


45 


-.0094 


-.0068 




-.2172 


-.2172 




-.2266 


-.2240 




50 


-.0275 


-.0200 - 


.0417 


-.1901 


-.1901 


-.1458 


-.2176 


-.2101 


-.1875 


55 


-.0409 


-.0300 




-.1301 


-.1301 




-.1710 


-.1601 




60 


-.0868 


-.0683 - 


.0509 


-.1156 


-.1253 


-.1944 


-.2024 


-.1933 


-.2453 


65N 


-.0868 


-.0868 




.0721 


.0721 




-.0147 


-.0147 




Wt. 

Avg 

(N. 

Hem) 


.0668 


.0581 


.0699 


-.1508 


-.1512 


-.1208 


-.0840 


-.0931 


-.0509 



Table VII. Zonally-averaged components R, -(E+Hp) and (Q^+S) of the 
surface heat budget. The averaged values as found in this study using 
subscript "l" for the 2/3-CL case and the subscript "2" for the full-CL 
case. Corresponding climatology values for the month of December re- 
ported after Berliand (1956) are denoted by subscript "3". Weighted 
averages are listed in the final line for the latitude range available 
in the N. Hemisphere. (All values in ly min *.) 
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204 - 




Figure 22. Distribution of R, - (E+H^ ) and of (Q v0 +S 0 ) 
as a function of latitude for the two cloud-model cases, 
(solid curves -- 2/3-CL case; dashed curves — full-CL case) 



130 



the 2/3-CL case, the corresponding is considerably larger 
in low- and mid- latitudes than R^, for the full-CL case. 

There is virtually no change in the - (E + H ) profiles over 
latitude (except poleward of 60N) . This means that the in- 
ternal heating rate of the ocean-mass, which is given by the 
sum of the two terms on the right side of (8-3), is generally 
greater in the 2/3-CL case than in the full-CL case. Here 
Q v0 is the required mean oceanic heat-flux divergence at 
latitude <j> . Figure 22 supports the conclusions of Von der 
Haar and Oort [1973] which requires increased lateral trans- 
port-divergences Q vq in the oceans, based upon the recently 
observed satellite -values of the global albedo. The values 

ascribed here to (Q + S ), seem to corroborate the satel- 

vx vo o' 1 

lite findings. 

The "hemispheric" mean value of (Q vq + S q )^ computed 
between 0 - 65N is 

(Q vq + S)^ = -.0839 ly min 1 

which corresponds to negligible cooling rate in the water-mass 
column, although it is five times as great as the correspond- 
ing atmospheric effect. 

C. OCEAN -TROPOSPHERIC HEAT -BUDGETS FOR ICE -COVERED OCEANS 
1 . Computation of (E + H p ) Over Ice-Covered Oceans 

The determination of E + Hp over ice-covered ocean 
was computed (following Equation A-36 of Kaitala, 1974)) 
with the following interface conditions used in Section V: 

E + H r = R + H. (5-11) 

T ice v ' 

H. = BBB (T -T ) . (5-10) 

ice v o g' 



131 



In latitudes where ice-covered oceans occur, R is negative in 
January, and H^ ce is usually positive but of an order of 
magnitude smaller than the magnitude of R. Hence the (E + 

Hp) values computed at both 4> = 60N and 4> = 65N from the 
170W meridian are both negative (here the Bering Sea is 
frozen on January 16, 1974). 

Solution for H^ ce and Hp + E at <f> = 60 and 65 N then gave 
the following results 





TABLE 


VIII. Surface heat- 


transports over 


ice at 170W 


Lat 




H. 

ice 


Ay min 1 


E + Hp 


Ay min 






2/3-CL 


FULL-CL 


2/3-CL 


FULL-CL 


4> = 


65N 


.0147 


.0147 


- .0721 


- .0721 


4> = 


60N 


.0041 


.0070 


-.0683 


- .0513 



The values of E + Hp of the foregoing display are used in the 
computation of zonal ly - averaged heat balance, both in the 
troposphere, and the ocean as reproduced in Tables VI and 
VII, respectively. 

The reason for the difference in computed values of E + H p 
associated with the 2/3-CL and full-CL cases at latitudes <J> 

> 60N is associated primarily with Equation 5-11 and the 
presence of ice. At <f> > 60N, the surface radiative quantities 
satisfy the inequality 

| R^ | < | R-p | , when CL > 0 

although both are negative. Furthermore, H. in either case 

1 C0 

is an order of magnitude smaller than R 2 or R^ , respectively. 
Hence, when CL > 0 
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|E ♦ H r | 2 < | E ♦ Hj, | j 



and there must be a smaller turbulent flux from air to ice 
in the full-CL case. The results of Table VIII at latitude 
60N exemplify this feature. 

2 . Budgetary Equations in Both Media with Ice-Covered 

Ocean 

Finally, H^ ce represents the result of a molecular 
conduction transfer rate, which leaves the ocean and enters 
the atmosphere over the ice by a process not included in the 
turbulent exchange rate E + . Hence at ice-covered grid- 

points the heat budget equations of the troposphere and 
water-body become: 

Air-mass over ice: 



The additive terms on the right side of 8-4, 8-5 did not 
modify the weighted-mean heating rates of the air-column or 
the water body significantly over the values resulting from 
Figs. 21 and 22. 




(8-4) 



Water-body under ice: 



Q + S = R - (E + H r ) - H . . 

wo o v T' ice 



(8-5) 
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IX. CONCLUSIONS 



The present work has been primarily involved with the 
development of the radiational two-layer cloud model adapted 
from Warner's earlier thesis [1974]. The data employed was 
oceanic data from the FNWC gridpoint analyses for 16 January 
1974. Here, it was found that the net radiative flux avail- 
able for the earth- troposphere system was best verified 
against January satellite climatology (Raschke et al. , 1973) 
when the Smagorinsky cloud-amounts were arbitrarily reduced 
by a factor of one-third. It is recognized that there were 
some minor differences in the comparisons which could be at- 
tributed to the one day space-time average as opposed to the 
average of several weeks for satellite data, and years inso- 
far as surface climatological data was concerned. 

Reduction of Smagorinsky cloud amounts by one third was 
found to provide more realistic global albedo estimates than 
resulted using unreduced cloud-amounts. In the comparison 
with satellite data, the radiative model with a reduced CL 
gave slightly higher albedo estimates in the tropics and 
mid- latitudes and slightly lower albedo in the high- latitudes . 
The radiative model still lacks an effective surface albedo 
model over ice. 

With the reduced cloud amounts greater surface net heating 
rates in the low to mid latitudes were computed than when 
using the original Smagorinsky cloud amounts. This latter 
finding is in accord with recent heat balance studies using 
satellite data which have found an increase in surface 
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heating particularly in the tropics. It is felt that further 
investigation of the tropical radiative balances are neces- 
sary. Thought should also be directed towards the geometry 
of solar radiation reflection problem among cumulus subgrid 
cells in the tropics as contrasted with the larger-scale two- 
cloud layer model used in this study, which proved fairly 
reliable in middle and high latitudes. Moreover, specifica- 
tion of the primarily cumuliform cloud distributions in the 
tropics is evidently closely related to the subgrid scale 
convection model rather than to the large-scale distribution 
of relative humidity (of Smagorinsky) . It is recommended that 
further testing of the radiative model for the different sea- 
sons of the year continue, particularly with the one- third 
or similar reduction in the Smagorinsky cloud-amounts. 

Further differences from climatological data might be 
attributed to the treatment of cloud reflectivity and absorp- 
tivity for solar wavelengths, as well as the parameterization 
of the sea-air exchanges of sensible and latent heat. In the 
latter case a new interpretation of the simple bulk formulas 
and of a redefined distribution of the turbulent sensible and 
latent fluxes in the planetary boundary layer parameterization 
model seemed to give a new, conceptually satisfying view of 
the planetary boundary layer model. Additional testing with 
the redefined planetary boundary layer or any new interpre- 
tation of such, might also prove productive in providing bet- 
ter comparisons with recent investigations such as the GATE 
and BOMEX experiments involving surface heat budget calculations. 
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APPENDIX A: EMISSIVITY e FORMULAS AS A FUNCTION OF WATER 

VAPOR AND CO p C ABSORBER MASSES ALONG PATH OF 
INTEGRATION. 



The general form of the calculation required for computing 
net IR flux is: 



where B nn = Stephen-Bolt zman blackbody flux at T, n and e 
is the emissivity of the atmosphere as a functipn of both 
water vapor and CC^ absorber masses along the path of 
integration. The hatched areas in Section III figures are 

®10 

represented by J e dB which is the key integral to be 

wc 



computed. In order to perform this calculation, Sasamori 
[1968] proposed emissivity formulas of the following form: 



for the layer 6-8 relative to the level 6 as reference level 
''or temperatures T ^ 210 deg k where £ v _ c is very nearly 
temperature independent. Also note that U is the water 
vapor absorber mass given by (2-7) measured from the 
reference level 10, 6, or 2 and C is the correspondingly 




(A-l) 



B=0 



e wc (6 ’ 8) = {,2it0 lo Sio [U(6 ’ 8)] + * 622} 

+ .07262((1. - .62556[U(6,8) + - 0286 ] * 26 



« [log lQ C(6 ,8) + 1.06^]} 



( A-2 ) 



136 



carbon dioxide absorber mass. Water vapor emissivity alone, 
contribution is represented by the quantity within the brace 
of Eq . A-2. The second brace indicates the additional contri- 
bution made by the part of the CO^ emissivity transmitted 
through the water vapor "overlap" near 15um, i.e., the part 
not already counted in this spectral region by the water 
vapor integration. 

Equation A-2 was used to compute atmospheric emissivities 

over all the layers (8, 10), (6, 10), (4, 10) (1, 10) 

for which the running mass values U(8, 10), C(8, 10), etc. 

were ready to be implied from the radiative soundings. To 

complete the F + integration given by (A-l) from level k = 1 

to the origin (B = 0, e =1.0), on the Yamato Radiation 

w c 

Chart (Appendix Pig. A-l), U(T^), C(T^) soundings were 
adopted as isothermal from T = T-^ (k = 1) to T = T q (k = 0), 
and then extended from the last "sounding" point (U[(0, 10), T-^]) 
along the last emissivity value for T < T^ = 210 deg k where 
the emissivity of the water vapor is known to be temperature- 
dependent. This value of emissivity is denoted e wc and is 
written, after Sasamori [1968] 

^[11(0,10),^] = 8.34 T.^* 353 log l° U(0}10) “ -W 

* U(0,10) C -* 03 ^5 logic U(0.10) - .705] 

* C 8 . 0/ (.353 log 10 U(0,10) + 3.56)] 

+ . 07262 { ( 1 . - .62556[U(0,10) + .0286] -26 ) 

* Clog 10 C(0,10) + 1.064]} ( A- 3) 
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Figure A-l. Depiction of e given by Eq. (A-3) on the 
Yamato Radiation Chart. Area above heavy dark curve repre- 
sents F I above level k=l from a radiative sounding, 
wc 
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Here the final term in the braces of Eq . A-3 is almost the 
same form of the CC^ correction to the total emissivity is 
in Eq. A-2. In Eq . A-3, however, this term refers to the 
full layer (0, 10), the first term is temperature dependent 
as well as dependent upon the summed absorber masses relative 
to the reference level k = 10. It can be seen clearly that 
Eq. (A-3) is usable with proper final parameters U(0, k), 

C(0, k), and T-j_ for the determination of the curvilinear 
portion of the downward flux contribution through level k 
for all of the reference level cases to be investigated below, 
that is k = 2, 6. 

It should also be noted that e wc has been integrated over 
the temperature T £ for which the water vapor emissivity 
curves are temperature dependent (see Appendix Fig. A-l) so 
that the result of Eq . A-3 is really the average of the 

4 

final curvilinear emissivity value multiplied by = St T-j . 
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APPENDIX B: FORMULAS USED IN DEPICTING THE DISPOSITION OF 

INCOMING SOLAR (F( A) ) FROM LEVEL k = 2 TO THE 
EARTH’S SURFACE 

These formulas are depicted by Figs. 7 , 8, and 9 in 
Section IV. They are used in this study to calculate amounts 
of insolation absorbed by the atmospheric layers and the 
earth's surface. For further explanation and key to these 
formulas consult Section IV. C. Dashed separation lines are 
used to divide those equations involving absorptions in 
successive layers e.g., (2,4), (4,6), (6,8), etc. 

1. Case (1,1) Overcast in Both High and Low Cloud Layers 
(see also Fig. 7) 



F4 4- = F(A) (1-.271[U(2,4) Sec Z]‘ 3 ° 3 
F4 + = F4l(RA(l)) 

A244- = F(A) .271[U(2,4) Sec Z]* 3 ° 3 
A24i = F4t.271[U(2,4) Sec Z] ,3 ° 3 



F6 = F4 4- ( l-RA ( l)-A ( 1 ) ) 

A 46 = F4 4- (A ( 1) ) (B-lb ) 

A2 6 = A241 + A24t + A46 

TD68 = l-.27l[U(6,8)5/3]' 3 ° 3 

303 (B_2) 

TD910 = 1-. 271[U(9, 10)5/33° J 



where TD is a term to denote the Manabe-Moller transmissivities 



for insolation beneath a cloud layer. In this event. Sec Z 



is changed to the mean slant-path sec z — 5/3 [after 
Katayama, 1966]. 



F84- = F6MTD68) 

F8+ = F84-(RA(2)) 

F6+ = F8+(TD68) 

F6 + 4- = F6t(RA(l) ) 

F84-4- = F6 4-4-(TD68) (B-3a) 

A684- = F6 4- - F84- 
A68+ = F8t - F6+ 

A684-4 = F6 4- 4- — F8 4-4- 
A68 = A684- + A68t + A684 + 



F94- = F8 4- ( 1-RA (2 ) - A ( 2 ) ) + F84-4 ( 1-A (2 ) ) 
A89 = F84-(A(2)) + F84-4(A(2)) 



(B-3b) 



F104- = F94-(TD910) 

FlOt = F10 4- ( a(G) ) 

F9t = FlOt (TD910) 

F94-4- = F9+(RA(2)) 

F104-4- = F94-4 (TD910 ) (B-3c) 

A9104- = F94- - F104- 
A910t = FI Of - F9+ 

A9104-4- = F94-4- — F1044- 

A910 = A9104- + A910t + A9104-4- 

A810 = A89 + A910 



2. Case (1,0) Upper Overcast Only (see also Fig. 8) 

A24l = F(A) .271[U(2,*0 Sec Z]‘ 3 ° 3 

A24t = F(A){l-.271[U(2,fJ) Sec Z ] * 3 ° 3 }RA( 1) 

?0 o (B-4a) 

x'{.271[U(2,4) Sec Z]' 3U3 } 

A24 = A24t + A24t 

F6+ = F(A){l-.271[U(2,4) Sec Z] * 3 ° 3 } [1-RA(1)-A(l) ] 

A46 = F(A)(l-.27l[U(2,i|) Sec Z]* 3 ° 3 }A(1) (B-4b) 

A2 6 = A24 + A46 



F101 = F64- (TD610 ) 

F10+ = F104- (a(G) ) 

F6+ = FlOt (TD610) 

F 64-* = F6+ (RA(1) ) (B-4c) 

F1011 = F6++(TD6lO) 

A610+ = F6+ - F104- 

A610+ = F10+ - F6t 

A6104-+ = F64-4- - F104-4- 

A6l0 = A610+ + A6104- + A610 + + 
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3. Case (0,1) Low Overcast Only (see also Fig. 9) 



A 26 + = F(A) ,271[U(2,6) Sec Z]* 3 ° 3 

A68l = F(A) .271{[U(2,8) Sec Z ] * 3 ° 3 -[U( 2 ,6) Sec Z]’ 3 ° 3 } 
F8+ = F(A) {1-.271[U(2 ,8) Sec Z]’ 3 ° 3 } 

F8+ = F8+[RA(2)] 

F6+ = F8t{l-.271[U(6,8) Sec Z]* 3 ° 3 } (B-5a) 

F2+ = F8Kl- . 271[U( 2 , 8) Sec Z]‘ 3 ° 3 } 

A68+ = F8t - F6+ 

A26l = F6+ - F2 + 

A26 » A26+ + A26+ 



F9 » F8+[1-RA(2) - A ( 2 ) ] 
A89 = F8+[A(2)] 



F10+ = F9+ (TD910 ) 

F10+ = F10+(a(G)) 

F9+ = FlOt (TD910 ) 

F9-K = F9 + ( RA( 2 ) ) 

F10 + 4- = F9U(TD910) 

A910+ = F9+ - F10+ 

A910+ = F10+ - F91 
A910-H- =• F94-4- - F10 + + 

A910 = A9104- + A910+ + A910 + + 
A610 = A68+ + A68t + A89 + A910 



(B-5b ) 



(B-5c) 
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SUBROUTINE FOR COMPUTATION OF SURFACE WINDS FROM SECTION II. 



u 



o 

CM 



— o 

ZK 



o 


— I s - 






OO 








00 o 


> 








h- o 








-J h- 






OO 


LU 






Li. • 


o o — 






oo in 


# o— 






> i 


>h- in — 






V • 


h-IZO 






Oh 


-j~< i cos: 






>o 


UJ 0 O UJ sQcQ 


— * 




► • 


o-tt in 






V— h- 


oo 






<< 


+ M-fM ID w 


5! 




-JO 


•— O II ♦# 


CD 




<<-« 


X + ll Z • 


X 




► v# >« 


h-LU O CM 


o 




h* • ► 


-JM-— O — 


— « 




ZGN 


ujvom— *i/)v 


X 




— '2! V 


O *Oh-MM 






oo <r 


^ O 1 2 *h- 


1— 


i 


'**" • r*H 


x\ iu—<a: _j 


-J 




UO^XrH 


h» ^ in oo m 


LLI 


— i 


OO *• • 


-JvO • 4* *0 


CSS 




2S • CM h* 


UJ CM CM 


'mm* 


• < 


*^inxu. 


uC^ * r— < • 




r-H v0 ' 


3: • — 


^LAr-'— O — 




>UL, 


1 O si* «~H 


^ — — IX UJ LL 


LL 




LU -JO ► fOX>hH t rO *\ 


\ 


OX' 


z: »hxiao i 


1 0C '^OZO 


O 


no ■ 


Mh •'in v0 *f r\j cm O rO 


o— 


OrH 4" 


h- < m— ' •£> u.jx >- oo ^ -vOQw 


iAw 


O •'w i 


O-J^.h-0 • II 


II II *00000000 


oo 


ini— 



0< <CO^X>-rM 4-CQ imQDO^Qw< 

— w II II II <C£ 

— >LLLUa II tl LU LLI LU II LL^LLO JOOl'JO 

c^mq:u-Oqoqou.^cq^>o>o^ll 



o 

o 



h- OO O 

-4 r\\LC\ ro 



liiM 



VG=V1*5. 1479165E-01 
VSPC= ( . 3*VG) + 2. 2 
RETURN 
END 



SUBROUTINE FOR COMPUTING TEMPERATURE AT LEVEL 'X* FROM EQ 
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SUBROUTINE FOR COMPUTING LATENT HEAT FLUX IN THIS STUDY FROM EQ 
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SUBROUTINE FOR COMPUTATION OF HEAT CONDUCTION INTO ICE AND SENSIBLE AND 
LATENT HEAT FLUXES OVER ICE FROM SECTION V.D. 
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